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Abstract 

* ^ ■ ^ phenomenological theory of the fluctuations of velocity occurring in a fully developed homogeneous and isotropic 

^»-^" turbulent flow is presented. The focus is made on the fluctuations of the spatial (Eulerian) and temporal (La- 

<~i , grangian) velocity increments. The universal nature of the intermittency phenomenon as observed in experimental 

measurements and numerical simulations is shown to be fully taken into account by the multi scale picture pro- 
posed by the multifractal formalism, and its extensions to the dissipative scales and to the Lagrangian framework. 
The article is devoted to the presentation of these arguments and to their comparisons against empirical data. 
In particular, explicit predictions of the statistics, such as probability density functions and high order moments, 
of the velocity gradients and acceleration are derived. In the Eulerian framework, at a given Reynolds number, 
they are shown to depend on a single parameter function called the singularity spectrum and to a universal con- 
stant governing the transition between the inertial and dissipative ranges. The Lagrangian singularity spectrum 
compares well with its Eulerian counterpart by a transformation based on incompressibility, homogeneity and 
isotropy and the remaining constant is shown to be difficult to estimate on empirical data. It is finally underlined 
the limitations of the increment to quantify accurately the singular nature of Lagrangian velocity. This is con- 
firmed using higher order increments unbiased by the presence of linear trends, as they are observed on velocity 
along a trajectory. 

Resume 

Une theorie phenomenologique des fluctuations de vitesse Eulerienne et Lagrangienne dans un 
ecoulement turbulent Nous presentons une theorie phenomenologique des fluctuations de vitesse dans un 
ecoulement turbulent pleinement developpe isotrope et homogene. Nous mettons I'accent sur les fluctuations 
des increments spatiaux (Eulerien) et temporels (Lagrangien) de vitesse. La nature universelle du phenomene 
d'intermittence observe sur les mesures experimentales et les simulations numeriques est completement pris en 
compte par les arguments developpes par le formalisme multifractal, et ses extensions aux echelles dissipatives 
et au cadre Lagrangien. Get article presente les predictions de cette description multifractale et les compare 
aux donnees empiriques. En particulier, des predictions explicites sont obtenues pour des grandeurs statistiques, 
comme les fonctions de densite de probabilite et les moments d'ordres superieurs, des gradients de vitesse et 
de I'acceleration. Dans le cadre Eulerien, a un nombre de Reynolds donne, nous montrons que ces predictions 
ne dependent que d'une fonction a parametres, appelee spectre de singularites, et d'une constante regissant la 
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transition entre les regimes inertiels et dissipatifs. Le spectre des singularites Lagrangien est relie a son homologue 
Eulerien par une transformation basee sur la nature incompressible, homogene et isotrope de I'ecoulement, alors 
que la constante restante est difBcile a estimer a partir des donnees. Nous montrons finalement que I'increment est 
inadapto a quantifier prcciscmcnt la nature singulicrc dc la vitosso Lagrangicnnc. Cola est confirmc par I'utilisation 
d'increments d'ordres superieurs non biaises par la presence de comportements lineaires, comme nous I'observons 
sur la Vitesse le long d'une trajectoire. Pour citer cet article : L. Chevillard et al., C. R. Physique volume (annee). 
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1. Introduction 



As a long-standing cliallcngc in classical physics [1,2,3,4,5,6,7,8,9], fully developed turbulence is an 
archetypical non linear and non local phenomenon. When a flow is stirred at a large scale L, typically 
the mesh size of a grid in a wind tunnel or the width of the blades in a von Karman washing machine, 
the input energy cascades towards the small scales without being dissipated according to the classical 
picture of Richardson and Kolmogorov [1,2]. One of the most important objectives in turbulence research 
is to understand the processes that lead to this very peculiar cascading structure of energy. As far as we 
know, recent theoretical progress made in this long lasting field comes from the systematic analysis and 
description of empirical (experimental and numerical) velocity data. 

The first experimental measurements of turbulent velocity were performed in the Eulerian framework 
and focusscd on the longitudinal velocity profile. In this context, the Taylor hypothesis allowed to interpret 
the time dependence of the measurements obtained in a wind tunnel behind a grid or in an air jet as a 
spatial dependence [7]. These experiments gave access to the longitudinal velocity increments, i.e. 

Mx)^(u(x + £)-u(x)).^ , (1) 
where u is the velocity vector and i a vector of norm and the respective structure functions are given 

by 

M„(^) = {\5tur) . (2) 
For such flows, the Reynolds number deflned as 

7ee = ^, (3) 



where cr — ^ {(Slu)^) is a characteristic velocity at the so-called integral length scale L and v the 
kinematic viscosity, can be considered as very large compared to unity. At these high Reynolds numbers, 
Kolmogorov showed, in a first seminal article [2] using a dimensional analysis, that the second order 
structure function M2(i) behaves as a power law, i.e. 

M2{i) = {{5euf) = a' = CK{ef/'i'/' , (4) 

where ck is the Kolmogorov constant of order unity [10,11] and (e) the averaged dissipation that will 
be properly defined latterly. Equivalently, the power spectrum, i.e. the Fourier Transform of the velocity 
autocorrelation function, follows the celebrated Kolmogorov law, 

Eik) oc CK(e)2/3fc-5/3 , (5) 

where the proportionality constant can be calculated [10,11]. As stated by Kolmogorov himself [2], the 
laws predicted (Eqs. 4 and 5) are only valid in a range of scales, called the inertial range tjk <^ i ^ L {ov 
equivalently ^ A: <C Vii') delimitated by the Kolmogorov length scale r/K ^ LTZe^^^ under which 
dissipative effects dominate the physics. These laws have been successfully compared to empirical data 
[7]. If velocity fluctuations were Gaussian, these predictions could be easily generalized to higher order 
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Figure 1. (a) PDFs of signed longitudinal velocity increments in Modanc data [17]. Represented scales (from top to bottom); 
\n{l/L) = -6.4137, -5.6028, -4.6645, -3.6411, -2.7501, -1.8598, -0.8685, 0.1226. All curves are arbitrarily vertically shifted for 
the sake of clarity and their variance is set to unity. The solid curves correspond to our theoretical predictions (see text 
in section 2.2.2 and Ref. [43]). (b) PDFs of Lagrajigian temporal increments from the ENS-Lyon experiment, for time lags 
(from bottom to top, symbols •) t/T = 0.07, 0.16, 0.35, 1 and from Cornell a<;celeration data (symbols o, from Ref. [25]). 
Also, the curves are displayed with an arbitrary vertical shift for clarity, the variance is set to unity at any scales, and the 
original axis for the acceleration PDF (o) has been shrunk by a factor 4. Solid lines correspond to theoretical predictions 
(see section 3.1.3 and Ref. [24]). 

structure functions. As observed experimentally, the probability density functions (PDFs) of velocity 
increments undergo a continuous shape deformation, starting from the integral length scale L at which 
statistics can be considered as Gaussian, down to the dissipative scales where the PDF is highly non 
Gaussian [12,13,14]. This phenomenon is a manifestation of the intermittent nature of turbulence, as first 
underlined by Kolmogorov and Obukhov [15,16]. We show in Fig. 1(a) the estimation of the PDFs of 
longitudinal velocity increments obtained in the giant wind tunnel of Modane [17] at a high Taylor-based 
Reynolds number TZx = 2500 (Eq. 26 provides a link between the large-scale Reynolds number TZe and 
TZ\). The curves are arbitraly shifted vertically for the sake of clarity. This continuous shape deformation 
of these PDFs motivated Gastaing et al. to build a statistical description of these longitudinal velocity 
fluctuations [12]. Moreover, by a simple visual inspection, we can remark that the PDFs are not symmetric. 
This is related to the Skewness phenomenon associated with the mean energy transfer of energy towards 
the small scales that takes place in the inertial range. We will come back to this point latterly. 

The first part of this article is devoted to review the predictions that can be made using both the so- 
called multifractal formalism [7] and the propagator approach [12]. This description, that will be shown to 
depend only on a parameter function X'(/i) and a universal constant TZ* (independent on the flow geometry 
and the Reynolds number), accurately reproduces the non Gaussian features formerly presented. 

More recently, experimental [18,19,20,21,22,23,24,25,26,27]) and numerical [21,22,24,28,29,30] data have 
revealed a similar phenomenon in the Lagrangian framework (see recent review articles [31,33,34]). La- 
grangian velocity is defined as the Eulerian velocity of a fluid particle at the position X(t), initially at 
the position X(<o) via the following identity 

v(X(to),t) =u(X(0,t) . (6) 

In this framework, the study of the Lagrangian velocity fluctuations focusses on the Lagrangian time 
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increment defined as 



5^v{t) = v{t + r) - v{t) 



(7) 



with V a component of the Lagrangian velocity vector v (Eq. 6). A similar dimensional analysis d la 
Kolmogorov would give a linear dependence of the second order structure function 



where T = L/a is the integral time scale and the respective Lagrangian Kolmogorov constant [40,41]. 
This corresponds to a Lagrangian power spectrum of the form 



Once again, these laws are valid only in the respective incrtial range r^^ <C r <C T (or T^^ <C w <C t^^). 
Unfortunately, these laws cannot be easily generalized to higher order structure functions because of the 
fundamental non Gaussian nature of the velocity fluctuations. We show in Fig. 1(b) the estimation of the 
experimental velocity increments PDFs at various scales obtained at ENS Lyon [21] and the acceleration 
PDF obtained at the university of Cornell [25] . Let us first remark that the acceleration can be seen as 
a Lagrangian velocity increment at a scale t much smaller than the Kolmogorov dissipative time scale. 
Again we observe a continuous shape deformation from Gaussian statistics at large scale, to long-tail 
acceleration statistics at vanishing scale. This is again a manifestation of the intermittency phenomenon. 

We will show in this article that a similar statistical description can be developed in a Lagrangian 
context. The free parameters are again the corresponding singularity spectrum and a constant, as in the 
Eulerian framework. The Lagrangian singularity spectrum will be shown consistent with the prediction 
derived from its Eulerian counterpart using a transformation, presented latterly. 

This article is devoted to the presentation of a phenomenological theory of the statistics of the Eulerian 
and Lagrangian velocity increments, from the integral length scale L (or from the integral time scale T = 
L/a), down to the far dissipative scales based on the so-called multifractal formalism [7]. In this approach, 
statistical properties in the inertial range are assumed and compared to empirical data. This includes the 
classical K41 predictions "fc-^/S" (-gq 5-) g^^^^ "w"^" (Eq. 9), but also the intermittent (or multifractal) 
corrections. From this descriptive analysis, we predict the statistical properties of the (Eulerian) velocity 
gradients and (Lagrangian) acceleration as functions of the Reynolds number and of the corresponding 
singularity spectra. In this sense, we can consider the multifractal formalism, properly generalized to 
the dissipative range, as a phenomenological theory, i.e. statistical properties are assumed in the inertial 
range (given some free parameters fully encoded in the so-called singularity spectrum), and are predicted 
in the dissipative range. This was already recognized in Ref. [35]. To this regard, multifractal formalism 
can be viewed as a standard model of turbulence, and, as far as we know, it is the only formalism able 
to reproduce accurately higher order statistics of velocity increments. The most important perspective 
would be to establish a link between this formalism (in particular the existence of a singularity spectrum) 
and the equations of motion (i.e. the Navier-Stokes equations). This is out of the scope of the present 
article. Let us also mention some alternative statistical formalisms that are able to reproduce the non- 
Gaussian nature of the underlying statistics, such as, among others, "superstatistics" [36], continuous-time 
random walk models [37], vortex filament calculations [38] and kinetic equations approach [39]. As an 
interesting application of this multifractal formalism, we reinterpret former measurements by the group 
of Tabeling [44,45], concerning the velocity gradients flatness as a function of the Reynolds number. We 
end with a general discussion of the skewness phenomenon and of its implications on the modeling of 
velocity increments in the inertial range. In the Lagrangian framework, we introduce higher order velocity 
increments designed to quantify accurately the singular nature of velocity and justify their use. 
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Figure 2. Behavior of tiie velocity increment flatness as a function of the scales in a normalized logarithmic representation. 
The universal constant entering in the normalization is set to TZ* = 52. Different symbols represent different Reynolds 
numbers: Tlx = 140 (DNS by E. Lcvcque [42]), TZx = 208,463,703,929 (Experimental Helium jet [46]), 7^A = 380 (Air 
jet measurements by C. Baudet and A. Naert), TZx = 2500 (Modane's wind tunnel [17]). A straight line of slope -0.1 is 
superimposed on the data in the inertial range. 

2. The Eulerian longitudinal velocity fluctuations 

2.1. Behavior of the flatness of velocity increments from experimental investigation 

We show in Fig. 2 the estimation of the flatness F = l^^g^^y^2 of the longitudinal velocity increments 
in various flow configurations and various Reynolds numbers. The first set of data has been obtained in 
a Helium jet [46] at four different Reynolds numbers: TZ\ = 208,463,703,929, where the Taylor-based 
Reynolds number TZx is proportional to the square-root of the large-scale Reynolds number TZ,. (sec; Eq. 
26). Also are displayed the results in an air jet experiment at the ENS Lyon at Kx = 380 of C. Baudet and 
A. Naert. The higher Reynolds number {TZ\ = 2500) comes from the wind tunnel experiment in Modane 
[17]. For comparisons are also reported the results coming from a classical direct numerical simulation 
(DNS) of E. Leveque [42] at a moderate Reynolds number TZx = 140 based on 256'^ grid points in a 
periodic domain. 

The flatness F as a function of the scale £ is represented in a logarithmic representation. Flatness is 
divided by 3, i.e. by the value of the flatness of a Gaussian random variable. The scales are renormalized 
by the length scale L estimated such that the power-laws that are observed in the inertial range are 
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indistinguishable. Then, ln(£/L) is itself normalized by the Reynolds number, i.e. ln(7?.e /??.*), where 
TZ* = 52 is a universal constant, that will be shown to be linked to the Kolmogorov constant (subsection 
2.1.2). A similar procedure is applied to ln(i^/3). 

In this representation, the Kolmogorov length scale can be seen directly as \n{r]K / L) / ln{TZe/TZ*) = 
—3/4. It delimits the inertial range, i.e. ln{i/ L)/ \n{TZe/TZ*) e [—3/4,0], and the dissipative range, i.e. 
ln(£/L)/ln(7^e/7^*) < -3/4. 



2.1.1. The inertial range 

Let us first focus on the inertial range. Over this range, we observe a universal (i.e. independent on both 
the Reynolds number and the flow geometry) power law {l/L)~^-^ of exponent 0.1. In the same range, 
the second order structure function behaves as a power law, i.e. {{SiuY) ^ (^/L)^/"^ (data not shown). 
Furthermore, we can see that the statistics at scales £ larger than the integral length scale L are consistent 
with a Gaussian process since the Flatness is very close to 3. Indeed the flatness is slightly smaller than 
3. This can be justified theoretically from PDF closures [47,48]. We will neglect this effect in the sequel. 
At this stage, for £ > L, a probabilistic description of the velocity increments is straightforward, i.e. 

SfU aS, where 5 is a zero average unit variance Gaussian noise and = {{5lu)^), independent on 

the scale £. Henceforth, the symbol stands for an equality in law, i.e. in probability. It means the the 

PDF and all the moments of the random variables on the two sides of the equality are equal. 

In the inertial range, we need a probabilistic formulation of this power law behavior. This is clear at this 
stage that a simple Gaussian modeling is not enough since it would predict a scale-independent flatness 
of constant value 3. Furthermore, the parameters that we will use must be universal since the observed 
power-law is universal. The main idea to build up a probabilistic description is to mix Gaussian variables 
and, for example, to use a Gaussian random variable with a fluctuating variance. This was proposed in the 
so-called propagator approach [12]. The form of the fluctuations of this stochastic variance will be given 
by the standard arguments of the multifractal formalism [7]. The corresponding non-Gaussian modeling 
consists in writing: 

law " 



5eu = S, (10) 

with h a fluctuating variable, independent on the unit-variance zero-average Gaussian noise S, and char- 
acterized by its distribution function 



In the limit of vanishing values of £, h and 'D^{h) gain the mathematical status of the Holder exponent 
and the singularity spectrum respectively. Let us stress that the /i-distribution 'P\^\h) (Eq. 11) is indeed 
normalized, and the range of integration [/iminl/imax] depends on the precise shape of the singularity 
spectrum. One of the main hypothesis of the multifractal formalism is to assume that (h) is independent 
on the scale £ [7]. In this case, assuming mmfill — T)^ (h)] = 0, a steepest-descent calculation (see Appendix 
A) shows that (recall that h and 6 are assumed independent) 

[ph+l--D'^{h)\ 

{\s,un=<^''{mj^ [j) ri'\h)dh^^^aP{ml^jj , (12) 

where {\S\p) ~ F (^^) j^lP-K and F is the Gamma function. Hence, the form of the density 'Pl^\h) (Eq. 
11) implies that the structure functions behave as power-laws, with a set of exponents ([^^u]*') ~ £'^^ 
linked to P^(/i) via a Legendre transform [7] 



7 



Cp = min [ph+1- V^{h)] 



(13) 



It gives in a straightforward manner the behavior of the flatness in the inertial range, i.e. F = 3(£/L)^*' . 
Let us first mention that in a K41 framework, without any intermittency corrections, the singularity 
spectrum is equal to 2?-'^^-'^ (1/3) = 1 and 'D^'^^{h) = -co for h ^ 1/3. In this case, using Eq. 13, wc can 
easily show that C^^"*^ = p/3, and because of the linearity of the (p function, the flatness is independent 
on the scale £. Clearly, the universal power-law behavior of the flatness in the inertial range as shown in 
Fig. 2, is the signature of the presence of intermittency in turbulence. 

In the literature, several models for V^{h) have been proposed. One of the most widely used is the log- 
normal approximation, giving the simplest quadratic form of the (h) spectrum including intermittency 
corrections, as proposed by Kolmogorov and Oboukhov [15,16]: 



2C2 



(14) 



In this case, the proposed stochastic modeling of the velocity increments (Eq. 10) has a simple probabilistic 
interpretation: with /imin = — oo and /imax = +oo, the velocity increments are modeled as a Gaussian 
noise 6 multiplied by a lognormal multiplicator a{i/L)^. This was the proposition made in the propagator 
approach [12]. It is easily seen that (p = cip — C2^, leading to F = 3{£/L)~ . In the sequel, wc will 
choose /imin = and ftmax = 1 because (i) a numerical integration of Eq. 12 shows no difference with the 
indefinite case, (ii) as we will see, the extension to the dissipative range implies /imin > — 1 (see section 
2.1.2) and (iii) rigorously, only Holder exponents greater than and smaller than 1 arc accessible when 
using increments. Experimental and numerical data as displayed in Fig. 2 show that the flatness behaves 
as a universal power-law of exponent —0.10 ± 0.01. This corresponds to C2 = 0.025, and C2 is called the 
intermittency coefficient. Also seen on empirical data (data not shown), (^3 w 1 [49]. Thus, in the sequel, 
we will take ci = 1/3 + 3c2/2 « 0.37, very close to its K41 prediction 1/3. This defines completely the 
quadratic singularity spectrum (Eq. 14). 

Another widely used singularity spectrum is the She-Leveque spectrum [50], 



P^(/i) 



-1 + 3 



1 + ln(ln(3/2)) 
ln(3/2) 



- 1 



{h-1/9) 



ln(3/2) 



(/i- l/9)ln(/i- 1/9) . 



(15) 



This spectrum is based on log-Poisson statistics [51] and yields C,p^ = p/9 + 2[1 — (2/3)^/'^]. This implies 
that F « 3{£/L)~'^'^^, a power-law behavior in good agreement with the empirical data shown in Fig. 2. 
More sophisticated methods based on wavelets [52,53] agree on the difliculty to discriminate between log- 
normal and log-Poisson statistics, and none of these approximations, for the moment, have been derived 
rigorously from first principles (i.e. the Navier-Stokes equations). For these reasons, we will use, in the 
sequel, the simplest lognormal approximation (Eq. 14) that compares well with empirical data, at least 
for low order moments (p with p < 6 and that gives a Gaussian distribution for the /i-exponent, that is 
easy to manipulate (see in particular in the stochastic modeling proposed in section 2.2.3). 

Finally, the proposed stochastic modeling (Eqs. 10 and 11) includes a functional form for the probability 
density fimctions of the velocity increments. A formal derivation of these PDFs, starting from the product 
of two independent random variables is reported in Appendix B. In simple words, we can see that the 
random variable Sgu (Eq. 10) is fully defined when is given the law of h (Eq. 11) and 6 (let us say 
Vs), that we assume, at this stage, independent. To simplify the derivation of the PDF VsguiS/u) of S^u 
(see Appendix B for a more rigorous derivation), let us consider only the absolute value of the velocity 
increment. We get 



\n\5e,u\ = lnc7 + /iln- +ln|(5|. 

1j 
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Given that the random variables h and In |^| are independent, the PDF of In \5eu\ is the convolution 
product of the PDFs of /iln |; and ln\S\, as it was noticed in Ref. [12], namely 

P,^^(lnM) = I p,^|,|(inM_/,inl)7?^j^^(/jlnl)d(/iln|) 



/ 



^inH(ln^-/ilny)7'^(/i)dft. 



Noticing that Pi„|5|(ln |^|) = \S\V\s\m) and Pi„^(ln^) = \Seu\V\s,u\{\Seu\), we finally get 



-h 



Vi'\h)dh , 



(16) 



which shows the expression of the PDF of the absolute value of the velocity increment. Without further 
assumptions, this derivation cannot be generalized to derive the PDF of the signed velocity increment, 
since we are taking at one point a logarithm, but it clarifies the fact that defining in probability a random 
variable (such as in Eq. 10) allows to derive the associated PDF. 

Actually, Eq. 16 is also true for the signed velocity increment, as shown in Appendix B, and we get for 
the velocity increment PDF Pseui^eu) at the scale £ the following form: 



/■/(max 1 / f\~^ 



Siu f £ 



Vi'\h)dh 



(17) 



where Vsix) = exp(— .x^/2)/\/27r is the PDF of a unit-variance zero- mean Gaussian variable. A numerical 
investigation (data not shown) of the PDF (Eq. 17) shows the characteristic continuous shape deformation 
associated to the intermittency phenomenon, in a similar way than observed in Fig. 1(a). Once again, this 
PDF is symmetric, i.e. Vstui^iu) = Vsiui—Sfu) and thus fails to describe the skewness phenomenon. Wc 
invite the reader to have a look at the sections 2.2.1 and 2.2.3 that provide a formalism able to reproduce 
the asymmetry of the PDFs, as observed in empirical data. We will see in the following section how to 
introduce the dissipative effects in order to obtain predictions on the statistics of the velocity gradients. 



2.1.2. The dissipative range 

The flatness of velocity increments, shown in Fig. 2, behaves in a very different way in the dissipative 
range, i.e. for scales ln{£/L)/ln{TZe/'R-*) < —3/4. First of all, no power-law is observed. Then, we can see 
a strong Reynolds-number dependence. This is very different from the behavior observed in the inertial 

range. In the limit of vanishing scales, the velocity increments can be Taylor expanded and we obtain 
6iu{x) = £dxu(x), i.e. a linear behavior as a function of the scale. This implies that the flatness tends to a 

Reynolds-number dependent function, independent on the scale £, i.e. F{£) -^^^jpp- The aim of this 

section is to understand and model the transition from the observed power-law behavior of the flatness 
in the inertial range, to this scale-independent behavior in the dissipative range. But flrst of all, let us 
have a look at the empirical data. 

In Fig. 2, we see that not all of the measurements exhibit at very small scales a scale-independent 
flatness. Only the DNS {Tlx = 140), the air-jet (7^A = 380) and the Helium-jet [Tlx = 208) have 
succedeed in resolving scales smaller than the Kolmogorov length scale. Indeed, the measurements are 
difficult because the Hot-wire size is usually of the order of the Kolmogorov length scale. If the Reynolds 
number is too high, dissipative scales are too small to be resolved (as for TZ\ = 463, 703 and 929). In 
the Modane's wind tunnel, the integral length scale is very large, and despite the large Reynolds number 
{TZx = 2500), rjK should have been accessible (which is of the order of the millimeter). Unfortunately, 
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the noisy local environment of the giant wind tunnel and the high temperatures reached (~ 60° Celsius) 
implied by the strong level of turbulence, prevented electronics from working properly at high frequency. 

Given the experimental limitations, we sec that the flatness underlies a rapid increase. This was the 
subject of Ref. [42] . The main underlying idea is the differential action of the viscosity. The multiplicator 
cr{l/L)^ fluctuates in the inertial range. In the dissipative range, up to a (random) coefficient, the mul- 
tiplicator becomes linear with respect to the scale in order to be consistent with the Taylor's expansion. 
The viscosity will regularize this inertial-range singular behavior at a scale that depends on the strength 
h of the singularity: the bigger is the multiplicator, corresponding to smaller /i-exponents, the smaller is 
the scale of regularization. Purely kinematic arguments, proposed in Ref. [42], show that the width of the 
so-called intermediate dissipative range [fy^,?7^] [35], where ry^ and 77+ are respectively the smallest and 
biggest dissipative scales, indeed depends weakly on the Reynolds number 

In (^Jp) - \/^e , (18) 

the Kolmogorov length scale rjK lying in between. In the representation chosen in Fig. 2, this prediction 
(Eq. 18) implies that the observed rapid increase occuring in the intermediate dissipative range should 
be steeper and steeper as the Reynolds number increases. This is what is qualitatively observed as long 
as experimental technics were able to reach these very small scales. 

This can be fully modeled in the context of the multifractal formalism. To do so, one has to come 
up with a dissipative scale that depends explicitly on the strength of the multiplicator, or in a more 
straightforward manner, on the exponent h. Paladin and Vulpiani [54] first proposed such a /i-dependent 
cut-off. Their reasoning was based on the local Reynolds number. For a scale lying in the inertial range, 
a fluctuating Reynolds number can be defined using a fluctuating characteristic velocity = a{£/L)^, 
leading to the local Reynolds number TZe = vgi/v. We can see that at the integral length scale, the 
associated Reynolds number T^-i, is unique, does not fluctuate, and is given by T^i, = TZg = ctL/v. This 
allows to define unambiguously a dissipative length scale: a scale 77 at which the local Reynolds number is 
of order unity TZ^ = TZ* = 0{1). The order one constant TZ* , is a free parameter of the formalism, a priori 
universal. It has been introduced phenomenologically in Refs. [46,55,56]. We will take it to be TZ* = 52 
and we will show that it is related to the Kolmogorov constant Ck- From there, one obtains directly the 
main result of Ref. [54] , namely 

1 

f TZ \ ''+1 

nih) = l(^^) . (19) 

Based on this result, Nelkin [57] showed the implication of a fluctuating dissipative scale (Eq. 19) on 
the modeling of the velocity increments for scales lying in the dissipative range. Within our approach 
summarized by Eq. 10, the velocity increment at a scale smaller than the dissipative scale r;(/i) will be 
modeled in the following way: 

law (20) 

where S, as in Eq. 10, is a zero-mean unit variance Gaussian random variable, cr^ = {(6lu)'^), independent 
of the random exponent h. This probabilistic modeling (Eq. 20) is consistent with the Taylor expansion 
of the velocity increment. It implies a model for the velocity gradient [57], namely: 

o law cr / r]{h)\''~^ 

9.n=-(^— j S. (21) 

Also we see that the proposed modeling for i > r]{h) (Eq. 10) and for £ < T]{h) (Eq. 20) is continuous at 
the dissipative length scale £ = r]{h). 
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As in Eqs. 10 and 11, to fully characterize in a probabilistic manner the random variable 6tu (Eq. 
20), we need to define the distribution of the exponent h. The proposition of Rcf. [57] is to take for the 
distribution of h at the dissipative scales I < r]{h) a scale-independent distribution (up to a normalizing 
function 2{£)) that also has to be continuous at the transition i = ri{h), namely 

where [h) is the same universal function entering in the /i-distribution of the inertial range (Eq. 11), 
that we will approximate to be quadratic (Eq. 14). 

As a general remark, as previously underlined in Ref. [35], the multifractal formalism becomes pre- 
dictive. Indeed, given the description of the inertial range (Eqs. 10 and 11), and in particular given the 
singularity spectrum 2?^(/i), we can predict the behavior of the velocity increments in the dissipative 
range (Eqs. 20 and 22) by a simple continuity argument. As a consequence, using Eqs. 19, 21, and 22, 
the 2p*^-order moment of the velocity gradients is given by following function of the Reynolds number: 

^ 2p(h-l) + l-I>g(h) 

= (9* ^/j; (5;) (23) 

where {5"^^) = -^f^, and 2(0) = //j (^) """"^ dh. The distribution of a Gaussian variable 5 is even, 
so all the odd moments of both velocity increments (Eqs. 10, 20) and gradients (Eq. 21) are vanishing at 
this stage. To include the skewness phenomenon, that requires further probabilistic modeling, we invite 
the reader to take a look at section 2.2.1. 

The explicit form of even-order moments of velocity gradients (Eq. (23)) allows us to derive two impor- 
tant predictions: the computation of the average dissipation and the dependence on the Reynolds number 
of the flatness of the velocity derivatives. The dissipation e is a key quantity in turbulence theory [7]. In 
an isotropic and homogeneous flow, the average dissipation is related to the second order moment of the 
velocity gradients as (e) = lby{(dxuY). Using Eq. (23), one obtains 

2(fe-l) + l-P^(h) 



'''•-''"''T?mL. '''' 

In the limit of large Reynolds number, a steepest-descent calculation (see Appendix A) shows that, up 
to an order one multiplicative constant, the former integral, is dominated by the term (TZe/TZ*)^'^ , where 
X2 = — niin/j ^ ^(^-1)+^-^ (^) ^ Because by construction, and for a wide class of singularity spectrum such 

that (via a Legendre transform, c.f. Eq. 13) Cs = 1? which is the case with the quadratic approximation 
(Eq. 14) with ci = 1/3-1- 3c2/2, we can show [7] that X2 = 1- This implies that (e) is independent on the 
Reynolds number. This verifies a basic hypothesis of Kolmogorov, namely the finiteness of (e) at infinite 
Reynolds number. A precise estimation of the integrals entering in Eq. 24 (see Appendix A) shows that 

(e) « ^'4- (25) 

As shown in Refs. [42,43], Eq. 25 gives a prediction of the Kolmogorov constant ck- Indeed, in the inertial 
range, neglecting intermittent corrections, the structure function is given by {{Sgu}^) = (7^{£/L)^^^ = 

/ .\2/3 

Ci<-((e)^)^/^. Using Eq. 25, it gives ck = [■Yg-j , showing indeed that this constant TZ* is related to 
the Kolmogorov constant Ck- Using TZ* = 52, one flnds ck = 2.33, which is slightly bigger than what 
has be measured on empirical data c/f « 2 [10,11]. One of the reasons of this apparent discrepancy is 
related to the fact that the constant TZ* is clearly linked to the presence of intermittency and the implied 
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extension of the intermediate dissipative range, whereas the definition of the ck assumes the absence of 
intermittency corrections. The prediction of the second order moment of the gradients allows also to link 
precisely the (large-scale) Reynolds number TZe to the Taylor based Reynolds number TZ\ in the following 
way: 

= . (26) 

Another important prediction is the dependence of the flatness of derivatives on the Reynolds number 
[57]. From Eq. 23, using a steepest-descent argument, we can show that, up to multiplicative constants 
of order unity (given in Appendix A): 

-777^ — 7^ ~ 3 — - With v„ =min ^ — 

Using the quadratic spectrum (Eq. 14), TV = 52, and Eq. 26 for the Taylor-based Reynolds number 

dependence, we obtain: 



(27) 



3 (^^^ ' = lAinl^^ = 0.937e^-36 , (28) 



which is consistent with empirical observations [58,59,60,61]. More precisely, we compare in Fig. 3 the 
multifractal prediction (Eq. 28) with various experimental measurements and numerical simulations as 
compiled in Ref. [61]. The present theoretical prediction reproduces quantitatively the empirical observa- 
tions. Let us stress that the quantitative dependence of the flatness of velocity derivatives on the Reynolds 
number (Eq. 28) is a genuine consequence of the intermittent nature of turbulence. In a K41 framework, 
this quantity would be Reynolds number independent. 

2.1.3. Full multi scale description 

In the preceding sections (2.1.1 and 2.1.2), we have seen how to model velocity fluctuations in respec- 
tively the inertial and far-dissipative ranges. It remains to write down a formalism, based on these two 
well known limiting ranges, that reproduces velocity statistics in the entire range of scales, including the 
intermediate dissipative range. 

A first naive idea to gather both the inertial and dissipative ranges in a unified description is to consider, 
at a given scale £ and Reynolds number TJ.,., the /i-cxponents. Within the range [/imin; ^max], if < h* , then 
the velocity increment lies in the inertial range; On the opposite, if /i > h* , then the velocity increment 
lies in the dissipative range. The transition occurs at the exponent h* defined by ^ = ri{h*), namely 

h*iW = -(l + ^^) . (29) 



\tl{1/L) 

Using both the laws of the velocity increments in the inertial (Eqs. 10 and 11) and dissipative (Eqs. 20 
and 22) ranges, we then obtain 

p(lt-l) + l-P(h) T>(h\ 



ra) 



dh+ - dh 



(30) 



where Z(t) normalizes the probability densities, namely 



"'^rLM ""-L (i) 

Unfortunately, this model does not compare well with ciinpirical data because the transition is too sharp 
(data not shown, sec [63]). Furthermore, the present modeling gives a continuous but not diflferentiable 
modeling of the velocity increment. 
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Figure 3. (a) Theoretical predictions of the velocity increment flatness obtained using Eqs. 32 and 33 with a quadratic 
singularity spectrum (sec text) at different Reynolds numbers, as proposed in Fig. 2. The higher is the Reynolds number, 
the steeper the increase in flatness that takes place in the intermediate dissipative range, (b) Direct comparison of the velocity 
increment flatness obtained at Modane TZx = 2500 (symbols) and from the model (solid line), (c) Comparison of various 
observations of the Reynolds number dependence of the flatness of the velocity derivatives with the present theoretical 
predictions (Eq. 28): F = 1.3611°^^^ [62] (Athmospheric measurements, dashed line), F = 0.9171°^^^ [60] (wind-tunnel 
measurements, dash-dotted line), F = 1.147^^■34 [61] (DNS, dotted line) and F = 0.9311°^^^ for the present prediction Eq. 
28 (solid line), (d) Comparison of various observations of the Reynolds number dependence of the skewness of the velocity 
derivatives with the present theoretical predictions (Eq. 36): S = —0.2&R,'^^^ [62] (Athmospheric measurements, dashed 
line), S = -0.337e° 0'5 [60] (wind-tunnel measurements, dash-dotted line), S = -0.327^°•" [61] (DNS, dotted line) and 
S = -0.1757e° "4 for the present prediction Eqs. 36 and 39 (solid line). 



A continuous and difiFerentiable transition, inspired by the interpolation function of Batchelor to model 
the second-order strncturc function [65], has been proposed in Rcf. [64] as the ciiitirc range of scales. In 
Ref . [43] , a slight modification of the proposition of Ref . [64] has been introduced in order to make the 
transition compatible with the far-dissipative predictions of Ref. [57] . It reads 
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ijj3(,5 with 



{if 



ir,{h)) 



(l-h)/2 ' 



(32) 



where the random variable 5 is again a zero-mean unit-variance Gaussian noise, a = {{Slu) ), and 



1 



{V{h)-l)/2 ■ 



(33) 



The normalizing constant is again such that J^^.'^J^ v'^\h)dh = 1, namely 



m = I 

Jh„ 



1 



r,{h) 



-I (D(/»)-l)/2 



dh . 



(34) 



Note that at a fixed given scale i, the velocity increment (Eq. 32) and the /i-distribution (Eq. 33) tend to 
the inertial description (given by Eqs. 10 and 11) in the limit of infinite Reynolds number. In the same 
manner, at a given finite Reynolds number, this proposed description tends to the dissipative predictions 
(Eqs. 20 and 22)) in the limit of vanishing scales £ — >■ 0. 

We show in Fig. 3(a) the theoretical predictions of the flatness as a function of the scales i for the 
different Reynolds numbers previously investigated in Fig. 2. To do so, we integrated numerically Eqs. 
32 and 33 using a quadratic singularity spectrum, with C2 = 0.025 and Ci = 1/3 + 3c2/2, /imin = and 
hmax = 1 and TZ* = 52. The Reynolds number TZe is obtained from the Taylor-based Reynolds number 
using Eq. 26. The proposed description based on Eqs. 32 and 33, reproduces the main characteristics 
shown in Fig. 2, namely, the universal power-law behavior in the inertial range and the rapid increase 
of the flatness in the dissipative range. Obviously, the theoretical predictions do not suffer from a lack 
of resolution and so with the proposed renormalization of the scales and of the flatness, all the curves 
tend to a universal plateau given by Eq. 28 when £ — >• 0. The rapid increase that takes place in the 
intermediate dissipative range is consistent with the kinematic prediction given in Eq. 18, namely the 
slope of this increase behaves as 1/ ^ln{7le/TZ*) in this representation. We can see indeed that the higher 
is the Reynolds number, the steepest is the increase. 

In Fig. 3(b), we compare more precisely the velocity increment flatness obtained in Modane's wind 
tunnel [17] to the multifractal prediction. The model does reproduce both the inertial and intermediate 
dissipative ranges, given the experimental limitations to reach the far dissipative range. 



2.1.4. Reinterpretation of the Tabeling's data as a non trivial effect of the dissipative physics 

As an example of the implications of the present theory, we reexamine in this section the observations 
of Tabeling and Willaime [45] . These authors investigated fully developed turbulence in a Von-Karman 
flow in gazeous helium at low temperature. Varying the pressure of the gas, they could span an unusually 
large range of Taylor scale based Reynolds numbers 200 < Rx < 2200. They measured the local velocity 
using a hot wire probe. Through the Taylor "frozen turbulence" hypothesis, they could access to the 
longitudinal derivative of the velocity, and its flatness F (Eq. 28). Up to Rx — 600, their results are in 
good agreement with previous literature and the present prediction (Eq. 28). Surprisingly, at Rx — 650, 
F presents a maximum, and goes down up to Rx — 1000, then raises again slowly (Fig. 4). Tabeling and 
Willaime interpreted their results as some evidence of a transition in turbulent flows. Some comments 
suggested that this behavior of F could be due to the finite size of the probe. However, as remarked by 
Tabeling and Willaime, this size (about lO/xm) is much smaller than the Kolmogorov dissipation length 
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Figijre 4. Flatness of the velocity derivatives as a function of the Taylor-based Reynolds number in a (a) linear and (b) 
logarithmic representation: (o) experimental data from a Von-Karman flow in gazeous helium at low temperature [45], (solid 
line) Multifrax;tal prediction for the velocity derivative flatness (Eq. 28), (stars) Multifractal prediction for the velocity 
increment flatness at a fixed scale lojL = 2.5 10~^. 

r] at B.\ = 650. Moreover, sueh a limitation was expected to yield a saturation of f at a constant value, 
not a well pronounced maximum. 

Indeed, following the present multifractal theory, identifying the "velocity derivative" with a finite 
difference at a constant length gives a maximum for the flatness as shown in Fig. 4. This is due to 
the rapid rise of the longitudinal velocity difference flatness in the intermediate dissipative range. This 
maximum occurs when the length ig coincides with the lowest scale of this intermediate range, which is 
much smaller than rj. Using Eqs. 32 and 33 with a quadratic singularity spectrum (Eq. 14), TZ* = 52, 
stars in Fig. (4) show the behavior predicted for F if £o/L = 2.5 10^"^, in reasonable agreement with the 
size of the sensor. 

However, the present theory cannot predict a further rise of F as observed experimentally. Also, the 
width of the predicted peak is much wider than observed. The present theory can explain some of the 
surprising features observed, not all. 

2.2. Consistent description of the skewness phenomenon 

2.2.1. General discussion on the skewness phenomenon 

This section is devoted to the modeling of the skewness of the velocity increments. As we just saw, 
modeling the velocity increment as a Gaussian random variable multiplied by a random amplitude (see Eq. 
32) cannot reproduce the asymmetric nature of the distribution of velocity increments since; the Gaussian 
random variable S, and its independence on the amplitude /3e lead to vanishing odd-order moments, i.e. 
Vp e N, {{6eu)'^P~^^) = 0. Nevertheless, keeping the same probabilistic description as in Eqs. 32 and 33 for 
the second order moment of velocity increments, allows us to predict in a consistent way the third-order 
moment of velocity increments if we use the Karman-Howarth-Kolmogorov equation: 

((M^) = -^(eK + 6.«^. (35) 

In Ref. [43] (see also Ref. [66]), we compared experimental data to the predictions obtained for {{Sf^u)^) 
using the Karman-Howarth-Kolmogorov equation 35 and the second order structure function ((J^u)^) 
obtained from Eqs. 32 and 33. Predictions and emprical data compares well in the whole range of scales 
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[43] without additional free parameters. In particular, the level of skewness in the inertial range is well 
reproduced and can be shown to be related to the universal constant TZ* . If we neglect dissipative effects 
in the exact relation Eq. 35, and the intermittent corrections on the second order structure function, it is 
easy to see that S{(.) = {{6eu)^) / is independent on the scale £ and can be further approximated 
to S* = —12/7^* = —0.23, in excellent agreement with empirical data [43]. Moreover, applying a Taylor 
development on both Eq. 35 and on the multifractal predictions (Eqs. 32 and 33) for the second order 
structure function, in the limit of vanishing scales, we get the following multifractal prediction for the 
third order moment of the derivatives: 



[2/1-1 -I?^(/i)] 



n* 



2(fa-2) + l--Z 
h+1 



"(h) 



(36) 



where is a negligible additive term, coming from the Taylor's development of the normalizing factor 
Z(£) (Eq. 34) [43]. This prediction of the third order moment of the velocity gradient (Eq. 36) using 
a Batchelor-Meneveau type of transition between the inertial and dissipative ranges (Eqs. 32 and 33), 
depends on the singularity spectrum , measured in the inertial range on empirical data, and on the 
universal constant TZ* . We show in Fig. 3(d) the numerical estimation of the multifractal prediction 
of the Skewness of derivatives (using Eqs. 32, 33 and 36) as a function of the Taylor-based Reynolds 
number, using the quadratic singularity spectrum (Eq. 14) and TZ* = 52. This prediction is compared to 
empirical data [60,61,62] as described in the figure caption and compiled in Rcf. [61]. We observe some 
dispersion between the three different empirical skewnesses, although the dependence on the Reynolds 
number seems to be universal. The difference in amplitude could be due to a lack of experimental and 
numerical resolution, and/or to a lack of statistical convergence. Thus, if the multifractal approach fails to 
predict the value of the skewness, it does reproduce accurately the Reynolds number dependence. Indeed, 
using Eq. 36, a steepest-descent calculation shows that the skewness of the derivatives behaves as a power 
law of the Reynolds number, i.e. 

ln(-5(0))/ln(7^e)~XS-l , (37) 



with 



Xs = mm 

h 



2{h-2) + l-V^{h) 



h + 1 



3 . 

- mm 
2 h 



2{h-l) + l-V^{h) 



h+l 



(38) 



Using a quadratic approximation for the parameter function (Eq. 14), we get Xs — ^ = 0.067. This 
power-law dependence on the Reynolds number has been already obtained by Nelkin [57] using a different, 
although related, approach based on the asymptotically exact relationship {{dxuY) = —2v{{d1u)'^). As 
shown is Fig. 3(d), further numerical estimations of relation (36), once rephrased in terms of Taylor-based 
Reynolds numbers using Eq. 26, leads to the following dependence of the velocity derivative skewness on 
TZy. 

^^^""^'^ - -0.1757^r''• (39) 



'-t-OO 



2.2.2. Modeling the velocity increments probability density function 

Indeed, consistent predictions for higher odd order structure functions are needed to predict the shape 
of the full velocity increment probability density function (PDF). Several propositions were made in the 
literature to account for the asymmetry of the PDF linked to the skewness phenomenon [12,43]. To do so, 
we must modify the noise 5 entering in the probabilistic formulation Eq. 32, and/or correlate this noise 
S with the exponent h. Asymmetric PDFs can be obtained if we change the Gaussian random variable 5 
to a non-Gaussian noise, still independent on the scale £ and on the multiplicator Pg. More precisely, it 
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was proposed in Ref. [12] to consider the random variable 6 as being a variable of density Vs{S) that now 
reads 

l+as- 



Ps{S) oc exp 



(40) 



2 V" ' "^VT+s^, 

where as ~ 0.18 is a universal constant, independent on both Reynolds number and scales. The main 
problem using this peculiar noise (Eq. 40) is that it leads to non zero average velocity increments. This 
could be fixed by introducing a scale dependent free parameter that centers the whole velocity increment 
PDF. Furthermore, to reproduce the non trivial behavior of the skewness in the dissipative range, we 
need to modiiy the parameter as, and to make it dependent on both scale and Reynolds number. In 
this spirit, still based on the hypothesis of independence of the two random variables h and 5, a general 
development of the PDF of 5 on a basis made of the successive derivatives of a Gaussian, called the 
Edgeworth's development, was proposed in Ref. [43]: 



+ 00 



1 r/" 



n=0 



(41) 



where the coefficients Xn{£) are functions of the scale £. As previously shown, the symmetric part (even 
terms) is well described by a Gaussian noise, which means that Ao(^) = 1 and X2n{^) = for n > 1. 
The coefficient Xi{i) is set to zero since, from Eq. 41, (Siu) = —a{j3()\i{tj = 0. Under these hypotheses, 
the third order moment is then given by {{5iuY) = —Qa'^{P'^)\^{t), which fully determines the coefficient 
\z{t) thanks to the Karman-Howarth-Kolmogorov equation (Eq. 35). Importantly, A3(^) does depend on 
scale and Reynolds number. As Eq. 35 is the only available constraint on A„, it is tempting (as a first 
approximation) to restrict the expansion to A3: \2n+i{^) = for n > 2. Additional statistical equations 
involving higher order odd moments of d^u would be needed to give the next A2„+i(^). This would require 
further modeling (primarily to get ride of pressure terms), which is outside the scope of the present 
work. Unfortunately, this crude approximation for the odd terms \2n+i leads to severe pathologies, such 
as negative probability for rare large events, and is not consistent with higher order statistics such as 
hyperskewness {{Stu)^) / {(Siu)'^)^^'^ (data not shown). To remedy for this weakness, Ref. [43] proposed 
to modify the variance of the Gaussian associated to the third term in the development (Eq. 41), in the 
following way: 

1 r _5V2_^^(^)^(^2_^)g-^V(2a^)" 



(42) 



where A3 is fully determined by the exact relation 35, and a an add- hoc free parameter, close to unity [43], 
aimed at describing higher order odd statistics. Then, from Eqs. 32, 33 and 42, the velocity increment 
PDF can be written, under the hypothesis of independence of d and h as (see Appendix B): 



'PstuiSeu) 



Seu 



(43) 



Then, using a quadratic singularity spectrum (Eq. 14), TZ* = 52 and the add-hoc coefficient a = VO.9 
(sec Ref. [43]), the predicted PDF (Eq. 42) successfully compares to empirical data [43] as shown in 
Fig. 1(a) for various scales. The shape of the experimental velocity increment PDFs, from the inertial, 
to the intermediate dissipative and far dissipative ranges is well captured by the present theoretical 
prediction (Eq. 42), consistently with the behaviors of the flatness and skewness of the velocity increments. 
Unfortunately, at this stage, it is not possible to motivate the choice of the additional free parameter a. 
To avoid having recourse to this parameter, we arc forced to abandon the hypothesis of independence 
of the singularity exponent h and the noise 5. A formalism that takes into account possible correlations 
between these two random variables is presented in the following section. 
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Figure 5. (a) Numerical estimation of tiic theoretical prediction of the velocity increment PDFs at various scales 
in the inertial range (Eq. 44). All PDFs arc of unit variance and the scales used are (from bottom to top) 
£/L = 10-1, 10-1-5, 10-2, 10-2-5, 10-3, 10-3-5, 10""', 10-''-5, 10-5. (b) Numerical estimation of the velocity increments 
flatness F{£) (o) and the skewness S{i) (□) from Eq. 44. Solid lines correspond to theoretical calulations (Eq. 53). (c) and 
(d) Pre-multiplied PDFs {Siu)'^Vsgu(Siu) and {Siu)^Vsgu{Siu), absissa are renormalized by {(S(u)'^) and ordinates axe 
renormalized such that the integrals are unity. 



2.2.3. Probabilistic modeling of the skewness in the inertial range 

It remains to give a consistent description of the asymmetric part of velocity increment PDFs in both 
the inertial and dissipative ranges without invoking an additional free parameter a in the PDF of the noise 
S (Eq. 42). In this section, we propose such a formalism for the inertial range, and leave the extension to 
the dissipative range to future investigations. The main idea is to correlate the exponent h and the noise 
5 in the multifractal description provided by Eq. 32. 

Let us first assume that the velocity increment can be written as a product of a Gaussian noise 6 and 
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an amplitude {i/L)^: S^u a{£/L)'^5, as in Eq. 32. The main difference with former assumptions is to 

let possible a correlation between h and S. The simplest way to deal with such a probabilistic formalism is 
to assume h and 5 jointly Gaussian. In this case, see Appendix B, the velocity increment PDF Vs^uiSeu) 
can be written as 



-h 



SfU ( I 



,h 



dh 



(44) 



where Vs,h{S, h) is the joint probability of the random variables 6 and h given by: 



cxp 



2(1 



{5-msY , {h-rUhY 2p(5 - ms){h - mh) 



(45) 

In Eq. 45, ms and <jg (rcsp. m/i and tr^) stand for the mean and variance of the random variable S (rcsp. 
h). The correlation coefficient p{£) = lies in the range [—1,1]. At large scale £ = L one has to recover 

0"^ = ((<5lw)^)- This implies ag = 1. 

Consistently with Eqs. 11 and 14, the multifractal formalism sets the mean and variance of the random 
variable h to: 

where C2 remains the intermittency coefficient (ca = 0.025 in empirical data, c.f. Fig. 2). Using the 
Legendre transform (Eq. 13), the /i-distribution leads to a set of structure function exponent 



Using Eq. 44, we can show that 

/oo 
xVs^u{x)dx = a 
-oo 



ms- p{£)^Jc2HL/£) 



The velocity increment statistics are of zero-mean if and only if 

ms = p{£Uc2HL/£) . 
For such parameters, we obtain the following velocity increment moments 

{5iu) = 



(47) 

(48) 
(49) 



[l + p2c2ln(i/£)] 



{{5(,uf) = -2a^ j pVc2ln(L/£) [Ap'c^ HL/£) + 3] . (50) 

The remaining free parameter is the correlation coefficient p(£) that will be fully determine by the Karman- 
Howarth-Kolmogorov equation (Eq. 35) and the Kolmogorov constant ck (or equivalently by 7^*) . Indeed, 
neglecting dissipative effects, Eq. 35 gives {{5i>uY) = — |(e)^. The average dissipation is given as (e) = 
aa^ /L, where a = 15/7?.* = 0.2885 (Eq. 25). Using the expression of {{5iuY) given in Eq. 50, we obtain 
for p{£), after solving a polynomial of third order, the following scale dependence 

A{a) 



with 



A{a) 



10 



50a + 25^25 + 4a2 



Vc2 ln(L/^) ' 

1/3 



2 [50q + 25V25 + Aa^\ 



1/3 



(51) 
(52) 
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Using the scale dependence of the correlation coefficient p{(.) (Eq. 51), we can see that the proposed 

description is consistent as long as the coefficient p{t) remains bounded in between -1 and 1. It corresponds 
to scales t/L < exp{- (a) / C2) = 0.94, using a = 15/TZ* and C2 = 0.025. We finaly obtain for the 
Skewness and Flatness factors 



We see that the proposed description (Eqs. 44 and 45), fully determined by the intermittent coefficient 
C2, the Karman-Howarth-Kolmogorov equation (Eq. 35, setting z/ = 0) and the Kolmogorov constant Ck 
(or equivalcntly 7?.*), gives consistent predictions (Eq. 53). Indeed, if intcrmittcncy is neglected, ^3 « |C2, 
the skewness is constant 5 ~ —0.23 and compares well with empirical skewness in the inertial range [43]. 
Let us also notice that the flatness does not tend to its Gaussian value 3 when I ^ L, but to 3.07. This 
could be corrected if a large-scale cut-off is introdiiccd in an ad-hoc way. 

In Fig. 5 are shown the results of the numerical integration of Eqs. 44 and 45, using C2 = 0.025 and 
a = 15/7?.*. In Fig. 5(a) we recover the characteristic shape deformation of the velocity increments PDFs. 
In Fig. 5(b), the velocity increments Skewness and Flatness are estimated and compared succesfully to 
the analytical calculations given in Eq. 53. We also display in Fig. 5(c) and (d) the second and fourth 
order pre-multiplied PDFs. We can see that the core of velocity increments PDF dominates in the second 
prc-multiplied PDFs, whereas the tails of the velocity increments PDFs contribute signiflcantly to the 
fourth-order one. 



3. The Lagrangian framework 

Recently, several experimental technics [18,19,20,21,22,23,24,25,26,32,27] and massive numerical com- 
putations [28,21,22,24,30,32,68] have been developed aiming following fluid particles along their trajectory 
in a fully developed turbulent flow. As stated in the introduction, a phenomenology similar to the Eulerian 
one can be developed in a Lagrangian context. The goal of this section is to present such a phenomenology 
introduced in Ref. [24] , that has been compared to a compilation of empirical data in [24,33] . 

3.1. Probabilistic formalism of the inertial and dissipative ranges 

3.1.1. The inertial range 

In the same spirit as in the section 2.1.1 devoted to the Eulerian framework, a probabilistic formulation 

of the Lagrangian velocity time fluctuations can be written down [24]. As observed experimentally, at 
the large integral time scale T, related to the integral length scale L = aT, the statistics of both velocity 
increments and velocity are very close to Gaussian statistics. Indeed, in stationary, incompressible, homo- 
geneous and isotropic turbulent flows, the Lagrangian velocity statistics can be related to their Eulerian 
counterparts [4] (see also Section 3.2.1). The statistics of the velocity time increments at a time scale r 
(Eq. 7) can be modeled once again as the product of two independent random variables: 

Srv''^a(^yS, (54) 

where (5 is a zero-mean unit variance Gaussian random variable and = {{Stv)"^) = {{Slu)"^) = "^{u"^}- 
The exponent h is independent on 6 and fluctuates according to the distribution law: 
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where, according to the multifractal formalism assumptions, the Lagrangian singularity spectrum 'D^{h) 
is universal (i.e. Reynolds number independent) and independent on the time scale r. As with Eqs. 10, 
11 and 13, this probabilistic modeling (Eqs. 54 and 55) is consistent with a power-law behavior of the 
Lagrangian structure functions (l^^iil^) ~ r^", the exponents are related to the Lagrangian singularity 
spectrum r>^(/i) via a Legendre transform: 

=min [ph + l- V^{h)\ . (56) 

Dimensional analysis (Eqs. 8 and 9) leads to a K41 description of Lagrangian turbulence neglecting 
intermittency. In this case, 'D^{\/2) = 1 and 'D^{h) = — oo for h 1/2. Experimental and numerical data 
actually revealed the presence of intermittency [21,24,28,30,33] that can be characterized by a quadratic 
singularity spectrum [24]: 

V^{h) = 1 - ^ , with = -+4 and = 0.085 . (57) 

2c2 2 

Note that the obtained value of ci is consistent with empirical data and with the dimensional prediction 
^2 = 1- As in the Eulerian framework, the bounds of the integration domain in Eq. 55 are /imin = 
and /imax = 1 (see Section 4.1). We will see in the following (Section 3.3) that the quadratic approxi- 
mation (Eq. 57) is too crude to describe the statistics of the acceleration and is not compatible with a 
quadratic Eulerian singularity spectrum via the Borgas' transformation (see section 3.2). Nevertheless, it 
gives a consistent description of the Lagrangian velocity statistics in the inertial range as well as in the 
intermediate and far-dissipative range [24]. 

3.1.2. The dissipative range 

As previously reported, the quadratic Eulerian spectrum (Eq. 14) has been extensively compared to 
a large set of empirical data, some of them displaying a large inertial range [tiktL] (in particular the 
Modane wind tunnel data, TZx = 2500). Sophisticated signal analysis procedures [53,69] concluded that 
this quadratic approximation cannot be distinguished from the data and this for both the increasing 
part (positive order structure functions) and the decreasing part (negative order structure functions) of 
the V^(h) spectrum. In the Lagrangian framework, such an analysis is much more difficult because the 
experimental technologies (silicon strip detectors [19], acoustic scatering [21], rapid cameras [26,27], etc.) 
are not as efficient as a well known hot-wire probe under the Taylor's hypothesis [7]. More fundamentaly, 
at a same Reynolds number, the width of the inertial range is much greater in the Eulerian framework 
than in the Lagrangian counterpart. Indeed, dimensional analysis predicts a Kolmogorov dissipative scale 
TriK proportional to LTZe and in turn (jTjj^/tjk ^ Ti-V^ S> 1. 

Let us now adapt the arguments justifying a fluctuating dissipative length scale in Ref. [54] to the 
dissipative time scales [79]. A local Reynolds number can be defined as TZr = v^t/u, where Vr is a 
characteristic fluctuating velocity at the time scale t. From Eq. 54, we obtain: 

T^r = (^) . (58) 

Reproducing the argumentation developed in Refs. [54,79], a h-dependent dissipative time scale can be 
deflned as the scale such that the local Reynolds number (Eq. 58) is of order unity: TZt^ ^ TZ^ ~ 0(1). 
A priori, the universal constant TZ^ is different from the corresponding Eulerian one TZ* . We thus obtain 

[79], 



21 



(59) 



If we neglect the intcrmittcncy corrections (in a K41 framework), then the exponent /i = 1/2 is unique, 
and we recover the Kolmogorov's dimensional prediction 



= T 



7^t 



(60) 



For a time scale r smaller than the dissipative time scale Tn{h), the velocity time increment can be Taylor 
expanded, i.e. 5rV = ra + o(r^), a being the acceleration. Like for the Eulerian velocity increments (Eq. 
20), we obtain the following stochastic modeling of the velocity time increments in the far dissipative 

range: 



(- law T 
OrV = a — 
r<rn(h) T 



Tr,{h) 



h-1 



(61) 



In analogy with the Eulerian framework (Eq. 22) [57,79], the distribution of the /i-exponcnts docs not 
depend on the scale r (up to the normalizing function Z{t)) and is a function of the Reynolds number: 



T<T„{h) Z{t) \TZ'^ 



l-P-^(h) 
2h+l 



(62) 



From Eqs. 59, 61 and 62, we can derive the expression of the high order moments of the acceleration, in the 
same way we derived the prediction of the moments of the velocity gradients in the Eulerian framework 
(Eq. 23): 

^ 2p(fe-l) + l-P^(fe) 



Z{0) 



Tie 



2)1+1 



dh 



(63) 



^"^"^ dh and (S^^) = Let us stress that the acceleration odd-order 



where Z(0) = J^- {^y^^ _ , _ 
moments are predicted to vanish: (a^''+^) = as observed in data. Finally, acceleration PDF can be 
expressed as a function of the singularity spectrum (h) and the constant TZ^ : 



h-1 i-p^(fe) 

2h + l 2h+l 



aT ( TZe 
TV 



h-1 
2h+l 



dh 



(64) 



where Vs{x) = e ^^/^/\/27r. A detailed discussion of the Reynolds number dependence of acceleration 
variance and flatness is provided in section 3.3. 



3.1.3. Full multi scale description 

As in section 2.1.3, we need an interpolation formula linking the velocity increments statistics in the 
inertial (Eqs. 54 and 55) and dissipative (Eqs. 61 and 62) ranges. An adapted Batchelor-Meneveau form 
inspired from the Eulerian framework (Eqs. 32 and 33) was proposed in Ref. [24]: 



5tV aPrS with /3x = 



where again the random variable J is a Gaussian noise of zero-mean and unit-variance, and 



(65) 
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with 




Z{t) = / 
Jt 



,-|(0^(h)-l)/7 • 



(67) 



1 + 



As in Eq. 33, the normahzing function Z{t) is determined by imposing 

Al^r 'Ph i^)'^^ = 1. The free 

parameter 7 entering in Eqs. 65 and 66 actually controls the transition from inertial to dissipativc physics: 
the bigger 7 is, the steeper the transition. In the sequel, we will use the Batchelor value 7 = 2, as in the 
Eulerian framework. Let us point out that we have found 7 < 2 when dealing with experimental signals 
and focussing on the second order log-cumulant [24] whereas the value 7 = 4 was used in Refs. [68,33] to 
describe the logarithmic local slope. Following Appendix B, the velocity increment PDF then reads (Eqs. 
65 and 66): 



We show in Fig. 1(b) the numerical estimation of the Lagrangian velocity increment PDF (Eq. 68) 
at different scales using a quadratic singularity spectrum (Eq. 57) and for two Reynolds numbers corre- 
sponding to the two sets of experiments {Tl\ = 740 for the acoustic scattering measurements at the ENS 
Lyon [21] and TZx = 690 [25] for the Cornell's sillicon strip detectors). We can see that the multifractal 
model predictions for the PDFs compared well to both sets of experimental velocity measurements [24] . 
We have used = 0.0753, cf = 1/2 + c^-, 7 = 1.08 and TZ'' = 30 to describe ENS Lyon PDFs. For the 
Cornell acceleration data, we have used Eq. 68 with C2 = 0.079, of = 1/2 + C2 , 7 = 1-3 and TZ^ = 30. 
As far as ENS Lyon data are considered, to obtain these free parameters, we have fitted the second order 
cumulant, i.e. the variance of In |(5rw|, over the whole range of scales (inertial and dissipative) , TZe, and 
T assumed known, and defining and 7 as the minimizers of the quadratic error of the theoretical and 
empirical second order cumulant [24] . We have further shown [24] that a similar fitting procedure on DNS 
Data {TZx = 140) leads to = 0.086 and 7 = 1.98, a value of the transition parameter 7 closer to its 
Eulerian counterpart (i.e. 7 = 2, Eqs. 32 and 33). The low value of the parameter 7 found in experiments 
can be interpreted as resulting from a low-pass filtering induced by the finitesize of the particule tracers. 

As for Cornell's data, we did not use the true acceleration PDF (Eq. 64) because the measured fiatness 
is much lower than (i) the one predicted by the present multifractal formalism and (ii) what is obtained in 
DNS (see the discussion in section 3.3). We chose in Ref. [24] to consider the Cornell's PDF shown in Fig. 
1(b) as a velocity increment PDF at a small time scale r = 0.0029r, corresponding to r = r,,^/8.62 in 
units of the Kolmogorov's dissipative time scale (Eq. 60). Then, the respective intermittency coefficient 
C2 and 7 are extracted from Cornell's acceleration data while minimizing the quadratic error of the 
theoretical and experimental PDF. 

In Fig. 6(a) is shown the behavior of the flatness of Lagrangian velocity time increments as a function 
of the scale r, for various Reynolds numbers (from R\ = 130 to TZ\ = 740) and flow configurations 
(DNS and experimental von Karman flows), as previously done in the Eulerian framework (Fig. 2). As 
expected for scales r greater than the integral time scale T, the statistics of velocity increments are close 
to Gaussianity. This can be deduced from the Gaussianity of Eulerian velocity increments for i> L, using 
an ergodicity argument (see section 3.2). 

In the inertial range, for time scales in the range —1/2 < ln(T/T)/ ln(7^e/7^''') < 0, an universal 
Reynolds-number independent behavior is observed. Recall that in the Lagrangian framework, at a given 
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Figure 6. (a) Flatness of Lagrangian velocity time increments as a function of the normalized scales t/T. Lower Reynolds 
numbers correspond to numerical data (7?^;^ = 130, 180, 290 from E. Lcvcque [23,24] and TZ\ = 400 from the Roma group 
[30,68]). The highest Reynolds number TZx = 740 has been achieved experimentally at the ENS Lyon [21]. Time scales arc 
renormalized by the integral time scale T and we use 7?.^ = 30. (b) Comparison of the proposed predictions (using Eqs. 
65 and 66), using a quadratic singularity spectrum (Eq. 57) with C2 = 0.085, ci = 1/2 + C2, 7 = 2 and TZ^ = 30 (solid 
line) to the numerical data at TZx = 290. (c) Logarithmic local slope of the Flatness of the data displayed in Fig. 6(a). 
(d) Theoretical logarithmic local slope of the Flatness obtained from a numerical integration of Eqs. 65 and 66, using a 
quadratic singularity spectrum (Eq. 57) with C2 = 0.085, ci = 1/2 + C2, 7 = 2 and Ti^ = 30, and for the various Reynolds 
numbers as shown in Fig. 6(a). We furthermore display (solid line, without symbols) a theoretical logarithmic local slope 
at a very high Reynolds number (TZe = 10^°). 



Reynolds number, the width of the inertial range T/r^^ ~ TZe is expected smaUer than its Eulerian 
counterpart L/r]K ^ TZl^^. Therefore, we may think that observing clear power-laws for Lagrangian 
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velocity fluctuations asks for higher Reynolds number than in the Eulerian framework. Indeed, for the 
range of Reynolds number available, no clear power laws are observed. This is confirmed when we display 
in Fig. 6(c) the local logarithmic slope where no plateau is observed in the inertial range at any Reynolds 
number. Thus, as far as Flatness is concerned, no clear power laws are observed. 

Scales ln(r/T)/ ln(7^e/^^) < ^1/2 correspond to the intermediate dissipative scales. Note that the 
extension of the intermediate dissipative range in the Lagrangian framework is much wider than in 
the Eulerian framework. This is related to the fact that, as we will sec, Lagrangian velocity is more 
intermittent than the Eulerian one, implying stronger fluctuations of the local dissipative time scale t^(/i) 
defined in Eq. 59. In this range of scales, we can also notice a Reynolds number dependence, that could be 
interpreted, as previously done in the Eulerian framework (Section 2.1.2), as a direct consequence of the 
fliictuating nature of the dissipative time scale. When the scale r tends to zero, wc observe a saturation 
of the velocity increments flatness to the corresponding acceleration flatness. This representation shows 
clearly the Reynolds number dependence of the acceleration flatness that will be further discussed in 
Section 3.3. 

As we see, estimating the intermittency coefficient is difficult since, contrary to the Eulerian framework, 
no strict power laws are observed. Nevertheless, using Eqs. 65 and 66, we can predict the behavior of 
Flatness over the entire range of scales, not only the inertial range. In Fig. 6(b), we compare the DNS 
data {TZ\ — 290) against the proposed formalism (Eqs. 65 and 66), using a quadratic singularity spectrum 
T>^{h), an intermittency coefficient = 0.085 (and cf = 1/2 + C2 ), the Batchelor value 7 = 2 for the 
transition and TZ'^ = 30. Thus, we find that Lagrangian turbulence is more intermittent than its Eulerian 
counterpart. The present formalism reproduces quantitatively the behavior of the flatness in the inertial 
and intermediate dissipative ranges. In the far-dissipative range; (when r ^ 0), this formalism seems to 
overpredict the value of the acceleration flatness (see the discussion in section 3.3). Given this limitation, 
the comparison between theory and empirical data is very satisfactory. 

We can see that, even if no clear power laws are observed, the present formalism gives a realistic picture 
of the Flatness at any scale. We can also see that, even if the multifractal formalism assumes the existence 
of power-laws in the asymptotic limit of very high Reynolds numbers, the predicted Flatness does not 
exhibit a clear power law at the finite Reynolds numbers under investigation. This can be clearly seen in 
Fig. 6(d) where the predicted logarithmic local slopes of Flatness, for the various Reynolds numbers given 
in Fig. 6(a), are shown. Theoretically speaking, the fact that no power laws are obtained in the model is 
mostly related to the; wide extension of the dissipative range implied by the strong level of intermittency 
that prevents from getting an extended inertial range. Another reason that explains why no clear inertial 
range power-laws are obtained in the model is that, for the largest scales of the inertial range for which 
the ratio t/T cannot be considered as small, the steepest-descent calculation (Eq. 56), that predicts 
power-laws, is not a good approximation of Eqs. 65 and 66. To this regard, studying the behavior of 
structure functions in a relative way [30,32,33], in the spirit of the extended self similarity, allows to 
weaken large-scale anisotropic effects and shows more clearly power-law behaviors. Nevertheless, when 
working at a very high Reynolds number, the model indeed exhibits a clear power-law for the Flatness. 
This can been seen in Fig. 6(d) where we superimpose (solid line with no symbols), as an illustration, the 
predicted logarithmic local slope of Flatness for TZ^ = 10^°. We see indeed the presence of a plateau in 
the inertial range. 

3.2. The Borgas' argument: linking Eulerian and Lagrangian intermittencies 

This section is devoted to establish a link between the Eulerian singularity spectrimi V^{h) and its 
Lagrangian counterpart 'D^{h) [24]. We will mainly recall the work of Borgas [79] and invite the reader 
to have a look at this reference for a detailed derivation. An alternative dimensional derivation of an 
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equivalent relationship has been also proposed in the literature [70]. Let us also mention more general 
arguments developed on kinematic bases [71,72]. 



3.2.1. Ergodicity principle 

Establishing a relationship between Eulerian and Lagrangian fluctuations requires as basic statement 
some principle of ergodicity. In simple words, we will admit that in an isotropic, homogeneous, incom- 
pressible and stationnary flow, the Eulerian average of a physical variable, obtained from summing up its 
realizations over space is equal to its Lagrangian average, obtained from summing up its values along the 
trajectory of a particle. This was first recognized and formalized by Tennekes and Lumley [4]. 

More formally, let us consider a one-point physical variable J" (e.g. velocity, dissipation, pressure and its 
derivatives, etc.) that does not depend on the scale. In the Eulerian framework, this variable depends on 
the spatial coordinates and time T^{x, y, z, t). In a Lagrangian description of the flow, this variable can 
be written as a function of the initial positions of the particles and time T^{xq, yo,zo,t). The assumption 
of incompressibility allows the following relationship 



V 



lim ^ [ [ [ T^{x,y,z,t)dxdydz = lim ^ [ [ [ T^{xo,yo, zo,t)dxodyodzo , (69) 

-^ + 00 V JJJy V^ + OO V JJJy 



that expresses the fact that an incompressible fiuid continues to fill the box as it moves around [4]. The 
next step requires the assumptions of isotropy and homogenity. When ensemble averaging the equality, 
this expectation can be taken inside the integrals by linearity. Then, homogeneity and isotropy allow us 
to remove these expectations from these integrals because of the independence over the space. Thus we 
are left with 

{J^^{x,y,z,t)) = (jr^(a;o, 2/0,^0,*)) • (70) 
The particular case of the observables (x.y, z.t) = cxp[ik.u(x,y, z,t)] and T'" {xQ,yQ, zo,t) = 
exp[ik.v{xo,yo, zo,t)] was treated in Ref. [4]. The expectation of these variables are the characteristic 
function the Eulerian u and Lagrangian v velocities. The equality of these characteristic functions im- 
plies the equality of the distribution and of all the moments of each velocity components, namely Vg, 
(m?) = {v'^). As observed in data, the Eulerian velocity has statistics close to Gaussian, which implies 
from the ergodicity principle the Gaussianity of the Lagrangian velocity. This has been checked in ex- 
periments as well as in simulations (see the Gaussian values of the Eulerian and Lagrangian velocity 
increments flatness shown in respectively Figs. 2 and 6). 

The idea of Borgas was to apply this ergodicity principle to the (scalar) observable dissipation e, i.e. 
T^{x,y,z,t) = exp[ike(x,y, z,t)] and T'"{xo,yo, ZQ,t) = cxp[ike{xo,yo, zo,t)], with as a main outcome 
the equality of the moments of the Eulerian dissipation and of the moments of the dissipation as seen 
by the particle along its trajectory. The next subsection is devoted to recall the Eulerian multifractal 
predictions for the moments of dissipation, and to extend these predictions to the Lagrangian framework 
[79]. 



3.2.2. Multifractal description of dissipation fluctuations 

For the sake of completeness, we repeat here the arguments developed in Ref. [79] . Multifractal predic- 
tions have been historically developed [7,15,16,76,77,78] for the coarse-grained dissipation over a ball of 
size £ centered on the position r: 

ee{r,t) = ^f e{r',t)dr' . (71) 

For a scale £ lying in the inertial range, a formalism similar to the one developed for the velocity increment 
can be written (Eqs. 10 and 11) using an exponent a and a singularity spectrum f^{a). Using the notations 
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of Ref. [79], we can write eg = {e){£/L)" ^ (equality in probability law), and the probability to get an 

exponent a at scale £ being given by 'Pa\a) ^ {(./ LY~^'^^°-\ When the scale t enters the dissipative 
range, we must take into account the fluctuating nature of the dissipative scale (Eq. 19) parametrized 
by the exponent a, i.e. i]{a) ^ Tie The choice of Borgas to apply the crgodicity principle to the 

dissipation was influenced by the fact that, in the limit of vanishing scales, the coarse-grained dissipation 
leads to predictions for the point wise moments of dissipation. We get, neglecting multiplicative (Reynolds 
number independent, but q dependent) constants. 



-3i 



-i)+i-/"(<.) 



,<;(a-l) + l-/^(a) 



lim(4) = (e«) ~ (e)' 7^e da ~ {e)'^TZe ' ' . (72) 

In a Lagrangian formulation, we will consider the average of dissipation during a time scale r: 

e,(ro,t) = - / e{r„,t')dt' . (73) 

At-t'\<r 

In a similar fashion as developed for the Eulerian description, the multifractal description of the statistical 
properties of Cr can be written down using the exponent k and the respective Lagrangian singularity 
spectrum /^(k) [79], i.e. Cr = (e)(r/T)'*~^ (equality in probability law) and the probability to get an 
exponent k at scale t being given by 'pIi\k) ^ {t/T)^^^ When taking into account the fluctuating 
nature of the dissipative time scale r^^(K) ~ TZe'^^^^'^'^\ we get 



r«max <,(^-l) + l-/^(^) 



lim(e?) = (e«) ~ (e)« / Ue dn {eYUe ^ (74) 

Identifying the leading order Reynolds number power laws exponents entering in the Eulerian and La- 
grangian predictions of the pointwise moments of dissipation (Eqs. 72 and 74), we get a relationship 
between the Eulerian f^{a) and Lagrangian /^(k) singularity spectra: 

n„)=-i«+(i+i»)/-(Ji^). (75, 

This relationship is the main result (Eq. 5.6) of Ref. [79]. 



3.2.3. Refined similarity hypotheses 

To investigate the implications of the Borgas' relationship between the singularity spectra /^'^ of 
dissipation and the singularity spectra V^'^ of velocity, we need to use a dictionnary [7] between the 
statistical properties of the coarse-grained dissipation and of the velocity increments. This is provided by 
the Refined Similarity Hypothesis (RSH) of Kolmogorov and Oboukhov [15,16]. In the Eulerian framework 
in an rf-dimensional Euclidian space, this hypothesis reads: 

{{5euy) ~ (ef o /i = I and ©^(/i) = /^(a) +d-l. (76) 

Similarly, in the Lagrangian framework, the RSH hypothesis reads 

{{Srvy) - (e«/2)T9/2 o /i = ^ and V^(h) = /^(k) + d - 1 . (77) 

Only recently the Lagrangian RSH hypothesis (Eq. 77) has been verified in numerical simulations [73,74]. 
Using Eqs. 76 and 77 to reinterpret Eq. 75, we easily derive the following relationship between the Eulerian 
V^ih) and Lagrangian 'D^{h) velocity singularity spectra (with d = 1) [24]: 

V^{h) = -h+{l + h)V^(-^ . (78) 
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Figure 7. Comparing Eulcrian and Lagrangian singularity spectra: Eulerian quadratic singularity spectrum (Eq. 14, thick 
solid line), Eulerian She-Lcvcque spectrum (Eq. 15, thick dashed line), the thin solid and dashed lines correspond to their 
Lagrangian counterpart using the Borgas' transformation (Eq. 78), dotted thin line correspond to the Lagrangian quadratic 
spectrum (Eq. 57). 

This relation Eq. 78 can be inverted 



To end this section, let us emphasize that the relation 78 is consistent with a non intermittent (K41) 
picture of turbulence. In this case, the Eulerian singularity spectrum is such that V^iX/i) = 1 and 
V^{h) = -oo if /i 7^ 1/3. Using relation 78, we obtain 2?^ (1/2) = 1 and V^(h) = -oo if h ^ 1/2. This is 
consistent with {(Srv)"^) oc r and more generally with Eq. 77. 

We represent in Fig. 7 various singularity spectra used in the present article, in both the Eulerian and 
Lagrangian framework, as it was done in Ref. [24,63] . First of all, we display the Eulerian quadratic (Eq. 
14) and the She-Leveque (Eq. 15) singularity spectra using respectively a thick solid and dashed lines. 
One can see that the increasing part of the spectra superimposes on each others. It means that, in a purely 
inertial description of turbulence, the quadratic and She-Leveque spectra are indistinguishable when using 
positive-order structure functions. The decreasing parts are distinct. More sophisticated signal analysis 
methods than computing velocity increment moments are required to investigate the decreasing part of the 
singularity spectra, such as the wavelet transform modulus maxima (WTMM) [52] , the inverse structure 
functions method [80] or the wavelet leaders [53]. The conclusion of these investigations is that, given 
the statistical limitations and linearization effects, experimental signals are consistent with a quadratic 
singularity spectrum (Eq. 14) for both the increasing and decreasing parts. 

We furthermore display on the same plot (Fig. 7) the three different singularity spectra we are using. 
This includes the quadratic Lagrangian spectrum (Eq. 57) represented using a thin dotted lines, and the 
transformed Eulerian spectra using the Borgas' relation (Eq. 78) displayed using thin solid and dashed 
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lines. One can see that the increasing parts coincide for the three cases, quantitative differences are 

shown concerning the decreasing part. Let us stress clearly that the Borgas' relation being nonlinear, 
a quadratic spectrum in the Eulerian spectrum is transformed into a non quadratic spectrum in the 
Lagrangian framework. This explains why the thin dotted and solid lines are distinct. In the sequel, 
we will see that the decreasing part has a strong influence on the statistics of dissipative quantities in 
Lagrangian turbulence such as the acceleration flatness (sec section 3.3), giving several arguments to 
discriminate these spectra when compared to empirical data. Another important remark can be made at 
this stage. The Borgas' transformation (Eq. 78), given the Eulerian spectra formerly introduced, predicts 
the existence of ^.-exponent greater than unity. This has strong implications on the singular nature of 
Lagrangian velocity fluctuations as quantified with velocity increments and the universal character of 
acceleration. We invite the reader to sections 4.1 and 4.2 for further discussions of this important point. 



3.3. Prediction of the variance and flatness of acceleration 



This section is devoted to the multifractal predictions of acceleration. We will mainly focus on the 
variance and on the flatness. The even order moments of acceleration are given by Eq. 63. Recall that 
(5^) = 1 by deflnition, we get for the acceleration variance: 
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Using a Gaussian approximation of the former integral (see Appendix A), we get 
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2h + l ' ^''^ 

and hp is the /i-exponent for which f^{h,p) is minimum (in the same spirit as it is presented for the 
Eulerian case in Appendix A). Once rephrased in terms of mean dissipation (using Eq. 25, we obtain 
cr^/T = {e)TZ* /15) and in terms of Taylor-based Reynolds number TZx (using Eq. 26), we obtain 

(a2)=ao(e)3/V-V2, (83) 

where uq is a remaining non-dimensional quantity that includes intermittent corrections, tabulated in 
various flow conditions in Refs. [20,75] 
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When using the quadratic Lagrangian singularity spectrum (Eq. 57) with TZ^ = 30, we get: 

ao = 0.64437e^-^^*^ . (85) 

Alternatively; when introducing the Eulerian quadratic spectrum (Eq. 14) into the Lagrangian frame 
using Eq. 78, and TZ^f = 6, we obtain: 

ao = IMQSnl^^^^ . (86) 
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Figure 8. Multifractal predictions for the variance and flatness of acceleration. Different curves correspond to different sets 
of parameters 'D^{h) and T^t; quadratic Lagrangian spectrum (Eq. 57) and 7?.t = 30 (dotted line), quadratic Eulerian 
spectrum (Eq. 14) with 7?,t = 6 (solid line), and the She-Leveque spectrum (Eq. 15) with 72.t = 2 (dashed line). In the 
last two cases, the Borgas' transformation has been used (Eq. 78). (a) Acceleration variance given by Eqs. 85, 86 and 87. 
(b) Acceleration flatness based on the numerical integration of Eq. 88. Symbols correspond to numerical data: TZ\ = 130, 
180 and 290 (A) and Tlx = 400 (V) (Fig. 6), symbols o and □ correspond respectively to the flatness of a<;celeration and 
pressure gradient from Rcf. [61]. 



A very similar result is obtained when plugging the Eulerian She-Leveque spectrum (Eq. 15) into the 
Lagrangian frame using Eq. 78, and 7^t = 2: 



ao = 1.76037^^ 



(87) 



These three predictions show similar dependence on the Reynolds number [79] , but the multiplicative pre 
factor depends strongly on the constant 'R) (Eq. 84) that is found itself strongly dependent on the shape 
of the singularity spectrum (see also the following discussion on the acceleration flatness). We show in Fig. 
8(a) the Reynolds number dependence of the factor ao (Eq. 84) for the three different sets of singularity 
spectrum 'D^{h) and constants 'R) (Eqs. 85, 86 and 87). 

Similar predictions can be derived for the acceleration flatness. This study will underline the limitations 
of the quadratic Lagrangian singularity spectrum (Eq. 57), to describe the fluctuations of the velocity 
increments. The acceleration flatness can be expressed as the following function of the Reynolds number 
(Eq. 63) 



T{a) 



(a2)2 



3Z(0) 



fftmax in^\ 



4(h-l) + l--p^(h) 

^''+^ dh 



2(lt-l) + l-P^(h) 

^''+' dh 



(88) 



The predictions of the acceleration flatnesses obtained by numerical integration of the integrals entering 
Eq. 88 are shown in Fig. 8(b). These predictions are compared to DNS data, the ones used in Fig. 6 and 
the ones provided in Rcf. [61]. Some quantitative differences are observed for both the Reynolds number 
dependence and the multiplicative pre factor that depends strongly on the value of Ti) . An analytical 
expression of the flatness (Eq. 88) can be obtained once again using a Gaussian approximation (c.f. 
Appendix A): 
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where 



■ mm 

h 



p{h-l) + l-V^{h) 



and bp 



cm 



(90) 



2/1+1 

and 6^{h,p) is defined in Eq. 82. The exponents Xp control the Reynolds number dependence of the 

acceleration flatness. They depend only of the shape of the singularity spectrum T>^{h). In contrast, 
the multiplicative pre factor depends on both the shape of T>^{h) and TZ'^ . When using the quadratic 
Lagrangian spectrum (Eq. 57) with TZ^ = 30, we get 

T{a) = 0.0115n\'^^ . (91) 

When introducing the Eulerian quadratic spectrum (Eq. 14) into the Lagrangian framework using Eq. 

78, and 7^t = 6, we obtain 

T{a) = 0.41077^^'^^ , (92) 
whereas for the Eulerian She-Leveque spectrum (Eq. 15) and TZ^ = 2, we get 

J"(a) = 1.11927e^-^i^ . (93) 

The three predictions (Eqs. 91, 92 and 93) we made using the three different sets of parameters iV^{h) 
and 7?.''') show different behaviors. First, we can see quantitative differences on the dependence on the 
Reynolds number. When compared to data, the quadratic Lagrangian singularity spectrum (Eq. 57) gives 
an exponent 1.73 much bigger than for the two other parameter sets. This seems to be consistent with 
data only at the lowest Reynolds numbers, but not at the highest Reynolds numbers where the predicted 
exponent leads to some overestimate of the flatness. On the opposite, the quadratic Eulerian singularity 
spectrum (Eq. 14) and the log-Poisson one (Eq. 15), once re-expressed in the Lagrangian framework using 
the Borgas' transformation (Eq. 78), give a Reynolds number dependence consistent with data, especially 
at high Reynolds numbers. Second, the multiplicative pre factors depend strongly on the the shape of the 
respective 2?^ and on Ti) . 

Let us stress that the analytical approximations given in Eqs. 91, 92 and 93 differ from a numerical 
estimation of the integrals entering in Eq. 88. In particular, when the quadratic Lagrangian spectrum is 
chosen, differences are quantitatively significant. This is clearly due to the range of integration [hmin \ ^max] 
that is finite in the numerical integration (i.e. /imin = and /imax = 1) whereas it is assumed infinite (i.e. 
hmin = —00 and /imax = +oo) in the Gaussian approximation in order to get simple analytical formula 
(thus avoiding any corrections given by the error function Erf) . Let us point out that if we instead use (i) 
the range /imin = and /imax = 1 is used when working with the Eulerian quadratic singularity spectrum 
further transformed using Eq. 75, and (ii) the range /imin = | and /imax = 1 when working with the 
She-Leveque spectrum (Eq. 15), then the analytical prc!di(;tions (Eqs. 92 and 93) are very close to the 
numerical estimations of Eq. 88. This underlines the limitations of the quadratic Lagrangian singularity 
spectrum (Eq. 57) to represent the statistics of acceleration. Indeed, additional analytical work shows 
that the exponent h that minimizes the function 9^{h, 4) (Eq. 82) and leads to xi (c-f- Eq. 90) is negative 
when using the quadratic Lagrangian singularity spectrum (Eq. 57). This would imply the existence 
of unphysical negative /i-exponents. This drawback does not exist when using the Eulerian quadratic 
(Eq. 14) or the She-Leveque spectra (Eq. 15). Once again, this tells \is that the quadratic Lagrangian 
singularity spectrum (Eq. 57) is not able to reproduce the acceleration statistics (see section 4.1 for more 
detailed discussion about the range of integration [hmin, /imax])- 



31 



4. Further discussions regarding the singular nature of velocity 

4.1. Comments on the integration bounds /imin and /imax- 

This section is devoted to present further discussions on the choice of the integration bounds [hmin, hmax] 
that has been made to define the probabihstic models of Eulerian (c.f. Eqs. 32 and 33) and Lagrangian 
(c.f. Eqs. 65 and 66) velocity increments. This requires a clear definition of the Holder exponents h. Let us 
stress that the dissipative cut-offs, as given by the fluctuating dissipativc length scale r]{h) (Eq. 19) and 
by the fluctuating dissipativc time scale T^(/i) (Eq. 59), imply that hmin > — 1 in the Eulerian frame and 
/imin > —1/2 in the Lagrangian counterpart. In the sequel, we will consider only inertial range fluctuations 
because; the mathematical theory of singularities is well established (sec Rcf. [52] and references therein). 
Our choice of the integration domain will be based on these considerations. In Ref. [68], in which only 
the She-Levcque spectrum is considered, other choices have been made. 

The Holder exponent h{xo) at a point xq of a singular signal u{x) (the Lagrangian velocity v{t) could 
be chosen without loss of generality) is defined as the biggest exponent h such that it exists a polynomial 
Pjv of order N and a positive constant C verifying 

\u{x) - Pn{x - xo)\ < C\x - xo\^'^''°'> . (94) 

Let us stress that Pm is the Taylor's development of order N of the signal u and typically. A'' < h{xQ) < 
+ 1. In other words, h{x[)) is the exponent that defines the first singular behavior entering in the Taylor 
development of u in the neighborhood of x^: 

u{x) = u{xo) + {x- xo)u(i)(xo) + ... + ■^^^^^^(^^(xo) + C'\x - xol''("°' , (95) 

where u^^^ = d^u/dx^ and C a constant. This shows that the velocity increment is only sensitive 
to exponents < ft. < 1 because it is orthogonal to constants only. Thus, we chose for the Eulerian 
framework hmin — and ftmax = 1 when dealing with the quadratic spectrum (Eq. 14) and ftmin = 1/9 
and /imax = 1 when dealing with the She-Leveque spectrum (Eq. 15). Similar choices have been made in 
the Lagrangian framework, namely /imin — and /imax = 1 when dealing with the Lagrangian quadratic 
spectrum (Eq. 57) and the Eulerian quadratic spectrum (Eq. 14) once transformed in the Lagrangian 
frame using Borgas' relation (Eq. 78), hmin = 1/8 and /imax = 1 when dealing with the She-Leveque 
spectrum (Eq. 15) once reexpressed in the Lagrangian frame. 

As far as the quadratic spectrum (Eq. 14) is considered, let us stress that the Eulerian predictions based 
on Eqs. 32 and 33 are similar if one uses /?,„iin = and /imax > 1 or /imin = ci — \/2c2 ~ 0.1472 and ftmax = 
Ci + \/2o2, ~ 0.5944, the two values of h such that {K) > 0. In the Lagrangian framework, considering 
the quadratic spectrum (Eq. 57), the predictions based on Eqs. 65 and 66 depend quantitatively on the 
range of integration. This is related to the fact that the decreasing part of the Lagrangian quadratic 
spectrum (Eq. 57) reaches the zero value at /i « 1. This would strongly suggest the possible existence of 
/i-exponents greater than unity, exponents that cannot be measured with an increment, as explained at 
the beginning of this section. This is thus very tempting to study Lagrangian velocity fluctuations with 
multi scale objects orthogonal to polynomials as it may appear in the Taylor development (Eq. 95) if 
singularities h are greater than unity. This is the subject of the following section in which higher order 
increments are considered, as they were in Refs. [24,63,81,82]. 

4.2. Wavelet analysis of Eulerian and Lagrangian fluctuations. 

In order to check whether it exists singularities h greater than unity, we perform a higher order in- 
crements study of velocity fluctuations. In the Eulerian framework, we will consider A''*^-order velocity 
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Figure 9. Velocity increment Flatness for various sets of empirical data and increment order. Different lines correspond 
to different velocity increments orders: first (solid), second (dashed-dot) , third (dashed) and fourth (dotted), (a) Eulerian 
velocity flatness using the data from Modane's wind tunnel {TZx = 2500). (b) Experimental Lagrangian velocity flatness 
TZx = 740 from Ref. [21]. (c) and (d) Numerical Lagrangian velocity flatness TZx = 290. 



increments defined as 

diu{x) = Sg^^'ulx) = u{x + () — u(x) 
VN>2, S^^Kix) = + £)- 6'-/^~^^u{x) . (96) 

With such defined S^^^u velocity increments, singularities < h < N become accessible (see the discussion 
provided in section 4.1. 
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We display in Fig. 9(a) the results of the statistical analysis of the first four velocity increments on 
the experimental signal obtained in the Modane's wind tunnel [17] that provides the highest Reynolds 
number TZx = 2500. We choose to display only the velocity increment flatness F^^\t) defined as 

FW{i) = ^^'^J^^"^'^ . (97) 

In the inertial range of scales, namely £ G [10~^'^L;i], one can see that the observed power-law is 
independent on the; increment order TV. This means that singularities h greater than unity are not observed 
on Eulerian fluctuations, as it could have been expected from the shape of the well accepted quadratic 
or She-Leveque singularity spectra. This has an important implication on the multifractal formalism: the 
Eulerian singularity spectrum V^ih) is indeed measurable with the first order N = 1 velocity increment. 

At smaller scales £ lying in the dissipative range, one can see a quantitative dependence on the increment 
order A'': as N increases, the transition from the inertial and dissipative ranges is steeper and steeper. This 
can be easily understood while generalizing the present probabilistic approach (Eqs. 32 and 33) to any 
iV^'^-ordcr velocity increments, with thus the same singularity spectrum X'^(/i). Indeed, the dependence 
on the order N comes from the behavior in the dissipative range of the iV^^-order velocity increment, 
which, to be consistent with the Taylor development of velocity, reads 

S<,''K{x) ^^J^'d^uix) . (98) 

The fluctuating nature of the dissipative scale rj{h) (Eq. 19) is independent on the way one is looking at 
velocity fluctuations and thus can be parametrized with the exponent h, although the constant 7?.* might 
depend on the order N. The probabilistic modeling of velocity increments, consistent with the Taylor 
development 98, follows: 

the random variable 5 being again a zero average unit variance Gaussian noise, (a^^^)^ = {{S^j^^u)'^) = 
(m^) (this can be easily obtained using the binomial theorem), and the /i-distribution being unchanged 

an given by Eq. 33, but, in this case, ft-min = and /imax = A'^. A numerical investigation of the A^^'^-order 
velocity increment flatness (Eq. 97) based on the proposed probabilistic modeling (Eqs. 99 and 33) shows 
the same type of behavior as observed on the Modane's wind tunnel velocity data and depicted in Fig. 
9(a), namely a A-independent power-law in the inertial range and an intermediate dissipative steeper 
and steeper as the order A^ increases (data not shown, see Ref. [63]). In the Eulerian framework, the flrst 
order increment allows to measure fully the singular nature of velocity, the A'^^-order increment appears 
as a tool that highlights unambiguously the dissipative range. 

A similar study can be performed on the Lagrangian velocity using the A/^^^-order velocity time incre- 
ments deflned in the following way: 

6rv{t) = 5i^\{t) = v{t + r) - v{t) 
VAT > 2 , Si-^^vit) = Si^-^^v{t + r) - Si^-^\{t) . (100) 

We display the statistical analysis of the experimental measurements (in Fig. 9(b), TZx = 740) and of 
the numerical simulations (in Fig. 9(c-d), TZ\ = 290) of the Lagrangian velocity using the four flrst 
JV^^-order velocity time increments (Eq. 100). The results and interpretations are very different from the 
Eulerian case (presented in Fig. 9(a)). We can indeed see that, in the inertial range (t € [10~^'^r;r] 



{N-h)/2 ' (^^) 
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for the experimental signal and r e [10~^''^T;T] for the numerical simulation), the -order flatness 

F*^^^(t) underlies a transition between the N = 1 and N = 2 velocity time increment. Recall that in 
the Lagrangian framework, no clear power- laws are observed (see the discussion of Fig. 6). We display 
nevertheless the slope (i.e. —0.34 = — ) that we would be observed if experimental data and simulations 
were performed at higher Reynolds numbers. As for = 1 velocity increments, higher order {N > 2) 
velocity increments do not exhibit clear power laws but still do exhibit a steeper behavior. We are thus 
led to the conclusion that high order velocity increments do behave differently depending on the order 
N: the level of intermittency for iV = 1 is lower than the one observed for N > 2. This can be seen 
as a proof of the existence of singularities h greater than unity but smaller than 2, since the inertial 
range is found independent of the order N as long as A'' > 2. It also says that the singular nature of 
Lagrangian velocity can be captured by an increment (5*^^^ of order A > 2, the first order increment giving 
a biased singularity spectrum, that underestimates the intermittency coefficient, because of the existence 
of singularities greater than unity. 

In the dissipative range, in a similar fashion than in the Eulerian case, the rapid increase of the flatness 
is steeper and steeper as the order A^ of the analysing increment increases. This can be obtained using a 
similar probabilistic modeling of the velocity increment, as it was generalized in the Eulerian frame (i.e. 
Eq. 99): 

law ^m^ws with = ^^^' iN-ny, ' (101) 

1 + (^) 

the expressions of the dissipative time scale r^(/i) (Eq. 59) and the /i-distribution (Eq. 66) being un- 
changed, although the constant TZ^ might depends on A^ and the range of integration being /imin = and 
h^^ = N, and (aW)2 = ((jf )^)2) = ((<5f )«)2) = ^{u'). 

We have thus shown that some singularities exponent h greater than unity and smaller than 2 exist 
in the Lagrangian frame, and not in the Eulerian frame. This observation is consistent with the Borgas' 
transformation (Eq. 78) that indeed does predict these smooth behaviors in the Lagrangian trajectories 
starting from the well accepted Eulerian singularity spectra (Eqs. 14 and 15). This has a direct implication 
on the Lagrangian acceleration. Indeed, it suggests that acceleration might not be the right quantity to 
look at from a multi scale perspective. More precisely, it appears to be difficult to understand properly 
acceleration statistics as the dissipative limit of velocity increments statistics in the inertial range as it 
can be done in the Eulerian framework. This is also tempting to interpret the non universal features 
observed in empirical acceleration as a manifestation of the biased behavior of a flrst order increment and 
underlines the importance of Lagrangian velocity higher order derivatives as unbiased quantities. 



5. Conclusions - Perspectives 

We have shown all along the article that a phenomenological theory of turbulent velocity fluctuations 
can be written down. Assuming the probability laws in the inertial range, we can derive explicit pre- 
dictions for the statistics of velocity gradients and acceleration that can be compared to experimental 
and numerical data. In the Eulerian framework, assuming a quadratic (Eq. 14) or a log-Poisson (Eq. 
15) singularity spectra, and a universal constant TZ* , we can predict accurately the statistics of velocity 
gradients, including the flatness (Eq. 28) and skewness (Eq. 39). 

In the Lagrangian framework, predictions, when compared to experimental and numerical data, are not 
as satisfactory as in the Eulerian framework. In particular, the remaing constant TZ^ is shown to depend 
strongly on the precise shape of the Lagrangian singularity spectrum 'D^{h). Furthermore, it is shown. 
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using the Borgas' transform (Eq. 78), that the two widely accepted Eulerian singularity spectra (Eqs. 

14 and 15) lead to quantitative differences when computing the acceleration flatness (Eqs. 92 and 93), 
whereas a quadratic approximation for ^^{h) (Eq. 57) leads to an irrealistic acceleration flatness (Eq. 
91) at high Reynolds number. Nevertheless, these significant differences appear as a way to discriminate 
in the Eulerian framework the two former singularity spectra (Eqs. 14 and 15). We also mention the 
intrinsic limitations of the (first-order) velocity increment to quantify accurately the singular nature of 
the Lagrangian velocity (section 4.2). We thus propose to study higher order velocity increments (Eqs. 96 
and 100) in order to quantify precisely the singular nature of the velocity fluctuations. We showed that 
indeed singularities greater than unity do exist in the Lagrangian framework. 
As perspectives, we can mention: 

To investigate further the miiltifractal predictions of the velocity derivatives skewness (see section 
2.2.1). The present formalism underestimates the skewness of derivatives when compared to data (sec 
Fig. 3(d)), although the predicted Reynolds number dependence is realistic and a lot scatter is found 
in experimental data. To that regard, the experimental resolution is a key issue. Specially designed 
numerical experiments could give some hints on the influence of the resolution on this quantity, as 
it has been done in Ref. [83]. 

- A theoretical explanation of the value of the universal constant TZ* (or similarly of the Kolmogorov 
constant ck) would be welcome. More precisely, we may wonder how this value depends on the shape 
of the singularity spectrum (h) . This issue is even more crucial in the Lagrangian framework since 
(see section 3.3) the constant TZ^ is very different if we work with different singularity spectra. 

- We have seen that the (first-order) Lagrangian velocity increment is not adapted to study the singular 
nature of the velocity along the trajectories. The acceleration appears to be thiis an ill-posed quantity 
and we could question its universality. Indeed, depending on the flow geometry (von Karman flow, 
wind tunnel, fully periodic DNS, etc.), a lot of scatter is found in data (see [75]). To this regard, the 
second-order time derivative of Lagrangian velocity cPv/dt"^ may be more independent on the large 
scale geometry of the flow than the first-order one. Performing such a measurement (or a numerical 
simulation) will be much more difficult because of the presence of noise. It would be very interesting 
to quantify precisely the Reynolds number dependence of the variance and flatness of this quantity, 
as it was initiated in Ref. [61]. 

- The derivation of the Eulerian singularity spectrum from flrst principles is still missing. Recently, 
some progrc;sses have been made while studying the dynamics of the velocity gradient tensor Aij = 
duj/dxi along Lagrangian trajectories [84,86,87,85,88,89,90,91,92]. In particular, in Refs. [90,91] is 
proposed a stationary process for the velocity gradient tensor that predicts in a realistic way the lon- 
gitudinal and transverse intermittencies. This theoretical investigation gives an Eulerian singularity 
spectrum consistent with the one observed on data (Eq. 14). 

- In this article, we focus on two-points quantities, namely the velocity increments, fully determined 
by the corresponding probability densities Vsiu{5e.u) and Vs^viSrv). Nothing is said on the long- 
range correlated nature of velocity. Indeed, it is shown in Refs. [22] and [25] that the acceleration 
amplitude is correlated over the integral time scale T. In a similar fashion, it has been known for a 
long time that dissipation is correlated over the integral length scale L [6,93]. Taking into account this 
peculiar correlation structure of velocity allows to propose stochastic processes (one dimensional in 
the Lagrangian case, three dimensional vectorial for the Eulerian case) able to mimic the behavior of 
velocity in the inertial range [22,94]. It remains to propose such stochastic processes able to reproduce 
velocity statistics in the intermediate and dissipative ranges and more generally, to give a stochastic 
representation of the differential action of viscosity as pointed out by the fluctuating nature of the 
dissipative scales (Eqs. 19 and 59). 

We acknowledge for fruitfuU discussions P. Abry, P. Borgnat, P. Flandrin. We thank Y. Gagne and the 
ONERA for providing the velocity signal of the wind tunnel of Modane, P. Tabeling and H. Willaime for 
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the experimental data shown in Fig. 4, C. Baudet and A. Naert for the air-jet Eulerian data {TZx = 380), 
L. Biferalc, A. Lanotte, F. Toschi for the Lagrangian numerical data {TZx = 400 shown in Fig. 6) and N. 
Mordant for Lagrangian experimental data {TZx = 740). 



Appendix A. Gaussian approximation and steepest-descent estimation 



The aim of this Appendix is to elaborate on a Gaussian approximation that can be made to compute 
analytically the integrals entering, for example in Eqs. 24 or 28. Let us consider, without loss of generality, 
the integral 
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with TZe +00. We first perform a Taylor's development of the exponent of the Reynolds number around 
the exponent hp that minimizes it, namely 
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h—hj. 



Then, the even order moment of velocity gradient (Eq. 23) can be obtained in the following approximate 
way (including the normalizing constant ^(0)): 
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(A.4) 



h=h2 



Moreover from Eqs. A.4 and 24, we get a precise estimation of the dissipation: 
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Using the quadratic singularity spectrum Eq. 14 with C2 = 0.025, we get 



- 0.97- — - , 
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(A.5) 



(A.6) 



in very good agreement with the approximation made in Eq. 25. Similar calculations can be made for the 
velocity derivative flatness, and we get (for the quadratic singularity spectrum Eq. 14 with C2 = 0.025): 



3 X 0.99 X 



0.18 



(A.7) 
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once again is very close agreement with the approximation made in Eq. 28 without taking into account 
these corrections. 



Appendix B. Velocity increment cind gradient probability density functions 

We follow here the presentation of a classical textbook on random variables [67]. Consider two random 
variables x and y, and a function g{x,y). We form the random variable z as a function of the two random 
variables x and y, namely 

z = 5(x,y). (B.l) 

Consider the domain Dz such that 

{(x,y)e£>,} = {5(x,y)<4 . (B.2) 
The distribution function Fz{z) defined as the probability P to have z < z is given by [67] 

Fz{z) = P{z <z} = P{(x, y) e D,} = jj^ r^Jx, y)dxdy , (B.3) 

where 'Px,y is the joint-density of the random variables x and y. Then, the probability density function 
Vziz) = dFz{z)/dz is given by the derivative, with respect to z, of the distribution function Fz{z). 

Consider now the random variable Siu formed as the product of two random variables written as 
Seu = afit{h)5, with the function IHtih) > definite positive (Eqs. 32, 65 or 44). In this case, the domain 
D^^u is easily obtained 

{(/3^(/l), 6) e Di,u} = {ft e [ftmin, /imax], 5 S] " 00, -^]} . (B.4) 

We get for the distribution function F^^^ of the velocity increments 

Fs.uiSeu) ^ f VsA^^h)dSdh . (B.5) 

The probability density function of the random variable 5e,u is readily obtained from Eq. B.5 taking a 
derivative with respect to 5tu, and we get 

which justifies Eq. 44. If furthermore, the random variables h and S are assumed independent, as in Eqs. 
32 and 65, then the joint density can be factorized, i.e. Ps,h{S,h) = 'Ps{5)Vh{h) and the distribution Eq. 
43 is justified. 
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